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ABSTRACT 

Sixteen  lakes  in  Mount  Rainier  National  Park  were  sampled  between  Jul} 
3  and  September  27,  1985.   All  of  the  sampled  lakes  were  analyzed  for  total 
dissolved,  labile,  and  nonlabile  aluminum;  major  anions;  major  cations; 
dissolved  silica;  dissolved  organic  carbon;  and  pH .   The  objectives  of  the 
study  were  to  determine  the  sensitivity  to  and  extent  of  acidification,  and 
to  provide  seasonal  variability  information  with  respect  to  dissolved 
aluminum  and  acid-base  chemistry  of  the  sampled  lakes. 

Mount  Rainier  lakes  are  typically  highly  sensitive  to  acidification 
although  any  existing  level  of  acidification  is  below  the  detection  limits 
of  the  available  models.   Lake  aluminum  concentrations  follow  theoretical 
predictions,  appear  to  be  controlled  by  an  aluminum  trihydroxide  solid 
phase  similar  to  natural  gibbsite,  and  do  not  indicate  acidification. 
Seasonal  variation  of  aluminum  is  a  result  of  variation  in  organically 
complexed  aluminum  concentration  responding  to  variations  in  dissolved 
organic  carbon. 
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1.0   INTRODUCTION  AND  OBJECTIVES 


The  term  "acid  precipitation"  has  been  used  to  describe  precipitation 
that  has  a  pH  below  5.6,  the  pH  of  distilled  water  in  equilibrium  with 
atmospheric  carbon  dioxide  (0.03°o  CO  at  sea  level)  and  that  has  elevated 
concentrations  of  sulfate,  nitrate,  or  chloride  (Henricksen,  1933).   Acidic 
precipitation  has  been  reported  as  a  serious  problem  in  the  Northeastern 
United  States,  Southeastern  Canada,  and  Scandanavia  (Driscoll,  1980:  Dillon 
el  al . ,  1030:  Henricksen,  1973).   The  increasing  aciditj  jf  precipitai  b  i 
in  Europe  and  North  America  has  largely  been  attributed  to  an  increase  in 
sulfuric  and  nitric  acid  aerosols  produced  by  metal  smelting,  industrial 
processes,  and  fossil  fuel  consumption.   Released  into  the  atmosphere, 
sulfur  and  nitrogen  oxides  react  with  water  and  oxygen  to  fo  :i    ,r      •  ; 
nitric  acid. 

For  an  aquatic  system,  the  effect  of  acidic  deposition  is  dependent  on 
the  capacity  of  the  soil  and  water  to  neutralize  the  acids  arid  resist 
changes  in  pH .  The  carbonate  system  provides  the  dominant  buffer  capacity 
for  water.   Aluminum  compounds  play  a  significant  role  in  soil 
neutralization.   Previous  work  has  shown  acid  precipitai  i  us  to  resulf  in 
increased  aluminum  concentrations  in  the  Northeastern  United  States  and 
Scandanavia  (Johnson  et  al . ,  1931). 

Areas  susceptible  to  acid  precipitation  generally  have  shallow  soils, 
are  underlain  by  unreactive  bedrock,  and  have  small  drainage  to  lake  surface 
area  ratios.   The  Oregon  and  Washington  Cascade  Mountain  lakes  are  classi- 
fied among  the  most  sensitive  in  the  United  States  (Omernik  and  Powe 
1932).   Due  to  their  sensitive  nature  Cascade  lakes  would  be  expei  I  '  ti 
respond  quickly  to  low  levels  of  acidificat  ion. 


Acid  precipitation  has  been  reported  for  the  Central  and  North  Cascades 
of  Washington  State  with  low  level  acidification  reported  for  dilute  lakes 
in  the  Alpine  Wilderness  Area  of  the  Central  Washington  Cascades  (Welch  et 
al.,  1935).   Preliminary  investigations  by  the  United  States  Geological 
Survey  (U.S.G.S.,  1986)  in  terms  of  potential  lake  acidification  in  Mount 
Rainier  National  Park  suggested  early  stages  of  acidification  in  some  lakes 
because  of  elevated  levels  of  aluminum. 

The  increased  acidity  of  Washington's  precipitation  appears  to  be  due 
to  an  increase  in  sulfuric  acid  aerosols.   The  primary  source  for 
anthropogenic  sulfur  in  Washington  State  is  the  Centralia  Power  Plant 
(59,000  tons/year)  (Logan  et  al . ,  1982).   The  ASARCO  smelter  near  Tacoma, 
which  closed  in  1935,  had  been  releasing  80.000  tons  year.   A  significant 
amount  of  the  acidity  in  rainfall  may  also  be  derived  from  .  '',-'"  •>''■''■ 
contributed  to  the  atmosphere  by  internal  combustion  engines  in  the 
urbanized  Southern  Puget  Sound  lowlands  around  Seattle  and  Tacoma  (Welch  md 
Chamberlain,  1981).   Prevailing  southwest  winds  carry  these  atmospheric 
constituents  into  the  Cascades.   With  respect  to  the  lakes  in  Mount  Rainier 
National  Park,  the  objectives  of  this  study  were  to: 

1)  Characterize  the  major  ion  and  aluminum  chemistry  bui.ldiny   i  ihi 

data  base  of  the  U.S.G.S.. 

2)  Provide  seasonal  variability  information  on  water  chemistry  data, 

3)  Determine  sensitivity  to,  or  extent  of  acidification,  with  respect 
to  major  ion  and  aluminum  chemistry. 


2.0   THEORETICAL  BACKGROUND 

The  acidification  process  may  be  described  as  a  titration  of  a 
bicarbonate  dominated  solution  with  a  strong  acid.   Because  of  their 
prevalence  in  anthropogenic  emissions,  sulfate  molar  concentrations  ([SO 

4 

2- 
])  or  sulfate  plus  nitrate  molar  concentrations  ([SO     -  [NO   ])  are  often 

4  3 

used  as  surrogates  for  the  strong  acid.   During  the  initial  stages  of  the 

acid  titration,  bicarbonate  formed  through  weathering  will  consume  the  acid 

and  thereby  buffer  the  system  against  changes  in  pH.   During  this  phase  '  t 

is  difficult  to  detect  acidification  from  changing  pH  because  of  the  maski  ig 

effect  of  the  bicarbonate  buffering  and  the  natural  pll  fluctuations  induced 

by  temporal  variations  of  p     .   Only  when  the  bicarbonate,  buffering  is 

2 
nearly  exhausted  does  pH  variation  due  to  acidification  occur. 

Alkalinity  is  a  capacity  measurement,  unlike  pH  which  is  an  intens i t\ 

measurement,  and  is   unaffected  by  variation  in  p    .   Having  greater 

2 
natural  stability  and  being  initially  more  responsive  than  pH ,  alkalinity 

provides  a  more  suitable  parameter  for  detecting  acidification.   The  amount 

of  acidification  is  relative  to  the  alkalinity  consumed  by  the  acid 

titration.   Acidification  therefore  is  estimated  by  the  difference  between 

preacidif ication  alkalinity  and  present  alkalinity.   An  accurate  measurement 

of  original  alkalinity  is  usually  not  available  and  must  be  empirically 

estimated  from  contemporary  water  quality  data. 

Chemical  acidification  models,  used  to  estimate  original  alkalinity, 

are  based  on  the  concept  of  geochemical  weathering.   Inherent  assumptions  in 

the  empirical  models  include:  1}  electroneutra! i fcy  or  ionic  balance  requires 

that  a  solution  must  have  a  neutral  charge;  2}  there  is  no  nonatmospheri 

(J  idigenous)   source  <>f  chloride  in  the  watershed;  3)  carbonic  acid 


weathering,  which  produces  equivalent  amounts  of  bicarbonate  alkalinity 
( [HC0o  ])  and  cations,  dominates  the  system;  and  4)  cation  concentrations 
remain  constant  at  low  levels  of  acidification. 

2.1   Electroneutral ity 


To  satisfy  the  law  of  el ectroneutrality ,  the  equivalent  concentration 
of  positively  charged  ions  must  equal  the  equivalent  concentration  of 
negatively  charged  ions  in  aqueous  solution.   The  electroneutrality  equation 
and  carbonate  equilibrium  equations  provide  the  foundation  for  theoretical 
and  empirically  derived  models  of  acidification.   The  suitability  of  a  model 
for  determining  acidification  is  dependent  on  the  applicability  of  the 
assumptions  made  in  the  theoretical  equations  and  empirical  simplifications. 
The  derivations  of  theoretical  and  empirical  acidification  models  are 
presented  in  Table  1. 


Table  1 


Derivation  of  Theoretical  and  Empirical 
Acidification  Models 


Electrqneutral ity : 

1.  [H+]+[NH4^]-2[Ca2i]+2[Mg2*]-[Na+]-[K^]  - 
[HC03~]+2[S042"]  +  [cr]  +  [N0~3]  +  [0Hj 

Assume:   No  [CI  ]  input  from  watershed  (sea  salt  correction),  [N0o  ]  and 
[NH   ]  concentrations  are  negligible,  and  pH  <  7 : 

2.  [H"]+2[Ca2"]+2[Mg2+]-[Na+]  +  [K"]  =  [HCO„  "] -2[S0  ,2  ] 

By  defining  E  [cations]  as  E  ([M.Jz.)  where  M,  represents  the  netal  cation 
and  z.  represents  the  ionic  charge: 

3.  [H+]+  E  [cations]  =  [HCO  ~]+2[S0  2~J 
Alkalinity  Definition : 

4.  [Alkalinity]  =  [HC03  "] -2 [CC^2  ]  +  [OH~]-[lT] 
Simplified  at  pH  <  7: 

5.  [Alkalinity]  =  [HC0o^]-[H~] 

Substituting  above  expression  (eq.  5)  into  the  simplified  charge  balance 
equation  ( eq  3) : 

G.    [Alkalinity]  =  E  [cations]-2[S0  2~] 

4 

For  predicting  original  alkalinity,  assume: 

9- 

Xutural  (Background)  [SO  "  ]  is  negligible. 

4 

7.  [Alkalinity]  =  E  [cations] 

Substituting  empirical  correlations  for  E  [cations]: 

8.  b  ([Ca2+]+[Mg2+])  =  E [cations],  and: 

9.  b'  ([Ca" "])  r  E [cations] ;   where  b  and  b'  are  constants.   Substituting 
above  equations  8  and  9  into  equation  7: 

10.  [Alkalinity]   b(  [r,\2]-  [Mg2'~]  ) 


11 .   [Alkalinity]  -  b' ( [Ca2~] ) 


Carbonate  _Equi 1 i  brium : 


12.  C02(g)  -HpO  =  *HoC00  ;  KR ,  where  K'H  is  Henry's  constant,  and: 

13.  H2C03  =  HC03  +H  ;  Ka^  where  Ka1  is  the  first  carbonic  acid 
dissociation  constant. 

Genesis  of  alkalinity  and  cations  from  congruent  dissolution : 

1-1.   MCO  (s)-H   =  M   ^HCO   ;  where  M  =  mineral  (metal  cation).   By  combining 

equations  12,  13  and  14: 
15.   MC03(s)+C02(g)  =  M2^-2HC03" 


* 


By  combining  equations  12  and  13,  and  defining  K   =  K  K  *  : 


16.    K1  =  [HC03"][H^]/p 


C°2   * 


Substituting  from  eq .  7  into   K   equation: 

17.  K1  =  I  [cations]  [H  ]  'l>CQ 
Taking  logarithms: 

18.  pH  -  log Z   [cations] -log  K  -log  pCQ 

2 
Or,  substituting  from  empirical  equation  9: 

19.  pH  =  log(b' [Ca2+])-log*Ki-iog  pco 
Letting   a  =  log  b'-log  K  -log  prn  : 

20.  pH  =  a  +  log  [Ca2* ] 


Sea  salt  corrections  of  lake  water  ionic  composition  are  based  on 
assumption  two,  that  no  indigenous  sources  of  chloride  exist  in  the 
watershed.   Thus  the  chloride  present  is  assumed  to  be  derived  from  sea 
salts.   It  is  further  assumed  that  the  ion  ratios  of  sea  salt  and  sea  spray 
relative  to  chloride  are  the  same.   To  maintain  electroneutral  i  t.y ,  both 
cation  and  anion  concentrations  are  corrected  for  sea  salt  contributions 
relative  to  chloride.   The  assumption  that  the  ratio  of  each  ion  in  sea 
salt  or  natural  wet  deposition  is  invariant  relative  to  chloride 
concentration,  however,  is  not  strictly  correct.   Both  potassium  and 
calcium  are  associated  with  terrestrial  dusts  (Gorham  and  Martin,  1984;  and 
Deither,  1979).   The  effect  of  terrestrial  dusts  is  usually  minor  compared 
to  total  ionic  concentrations  (Wright  and  Gjessing,  1976). 

2.2   Empirical  Models 

Surface  water  acidification  is  estimated  by  comparison  of 
measured  data  with  the  theoretical  and  empirical  models  presented  in 
Table  1.   Equations  1-7  are  rigorous  if  the  assumptions  inherent  in 
their  derivations  are  verified.   Equations  8-12  involve  empirical 
simplification  but  may  be  applied  over  limited  conditions  where 
relative  ratios  of  cations  and  anions  do  not  vary  drastically. 
Alternatively,  substitution  of  carbonate  equilibria  (weak  acid  and 
solid  dissociation  reaction)  into  the  electroneutral ity  -  based 
equations  leads  to  pH  vs.  cation  concentration  equations  (e.g.  17- 
20,  Table  1) . 

Graphical  presentation  of  water  chemistry  data  facilitates  comparison 
with  theoretical  or  empirical  models  and  can  be  used  to  estimate 
acidification.   The  electroneutral ity  equation  applied  to  lakes  dominated 


by  carbonic  acid  weathering  results  in  a  linear  equivalence  between 
bicarbonate  and  total  cation  concentrations ,  i.e.,  slope  -  1  and  y 
intercept  =  0  ( eq  .  7,  Table  1).  A  slope  lower  than  unity  indicates  input 
of  cations  from  mechanisms  other  than  carbonic  acid  weathering  or  loss  of 
alkalinity  and  is  interpreted  as  evidence  of  acidification.  The  error  in 
the  y-intercept  term  is  usually  negative,  and  may  be  used  to  indicate  the 
concentration  of  ions  assumed  negligible  in  the  model  formulation  (Kramer 
and  Tessier ,  1982) . 

Empirical  simplification  of  the  el ectroneutra 1 i t y  equation  leads  to: 
[Alkalini ty]=b( [Ca  ]+[Mg2H ] )  and  [Alkalinity]=b' [Ca  ]  equal  ions  (Table 
1).   The  multiplier  (b  and  b')  terms  are  greater  than  unity  Lo  account  foi 
the  contribution  of  neglected  cations.   Variability  in  b  or  b'  is 
introduced  by  differences  in  cation  concentrations  betwpen  watersheds  and 
by  possible  acidification. 

A  rigorous  approach  for  constructing  a  pH  vs.  cation  diagram  would  be 

to  use  equation  18  (Table  1)  to  represent  carbonic  acid  weathering  since 

this  equation  is  theoretically  derived.   To  reduce  error  associated  with 

Prn  variation  the  theoretical  line  representing  carbonic  acid  weathering 
LU2 
is  drawn  for  a  fixed  p    and  the  pH  measurement  of  the  water  to  b^   plotted 

should  be  measured  at  that  p    or  be  corrected  for  variation  in  pCO„ 

2 
using   theoretical  equilibrium  equations.   Any  point  falling  above  the 

theoretical  equilibrium  line  indicates  that  sources  of  hydrogen  ions 
other  than  carbonic  acid  are  contributing  to  the  weathering  process. 
Henricksen,  (1980)  and  Kramer  and  Tessier  (1984)  applied  the  empiri- 
cally simplified  charge  balance  equation  between  alkalinity  and 

calcium  (eq.  11,  Table  1)  to  obtain  the  pH  vs.  [Ca"  ]  plot  using  the 

o  i 
equation,  pH  -a-  log  [Ca   ]  (eq.  20  Table  1).  High  variability  is 


introduced  into  the  constant  a  by  its  component  parameters  b 


:o 


a  lesser  extent  by  the  composite  dissociation  constant  *K  r>>:  carbon 
acid  (Table  1).   Tin-  empirical  constant,  h,  varies  sig:!:f'  i.nt]\  bel  eon 
surface  waters.  The  p  _   may  vary  by  one  order  of  magnitude  from  one 

9 
Ct 

location  to  another  or  at  the  same  location  with  time  (Stumm  and  Morgan . 
1981).   *K  is  affected  by  both  temperature  and  ionic  strength.   Despite 
the  simplifying  assumptions,  pH  vs.  ion  concentration  plots  are  often 
employed  to  determine  the  relative  acidification  among  lakes  since  many 
surface  waters  can  be  positioned  on  the  same  graph  for  comparison. 
Assumption  four,  that  cation  concentrations  remain  a  able  tvitl 
increasing  acidification,  has  been  questioned.  Low  level  acidic 
deposition  increases  weathering  reaction  rates  and  can  leach  cations 
rather  than  titrating  bicarbonate  (Galloway  et  al.,  1983).   This  obscures 
the  titration  ef  f  ec1  and  ;  s  u   eval  en1  in  grunridwi  '    '  ii  i  ted 
(Henricksen,  1980).   The  titration  concept  appears  most  directly 
applicable  to  watersheds  overlying  unreactive  bedrock  characterized  b\ 
shallow  underflows  and  short  retention  times. 


2.3  Aluminum  Chemistry 


When  primary  buffering  systems  are  exhausted  and  pH  declines  from 
between  6-7  to  below  5,  the  buffering  system  consisting  of  aluminum 
hydrolysis  compounds  becomes  dominant  (Driscoll  et  al . ,   1980).   At  pH 
values  below  5,  mineral  dissolution  reactions  increase  free  and  total 
i]  um.inum  cone  en  trat  i  ons  . 

Sensitive  watersheds  are  generally  underlaid  by  bedrock  that  has 
small  pools  of  readily  exchangeable  (weat!  Liable)  basic  ca  ions  and 


inability  to  retain  acidic  anions,  such  as  SO    .   Therefore,  if  high 
loadings  of  acid  are  not  attenuated  in  the  soil,  but  are  instead  removed 
by  drainage  water,  the  anionic  solute  must  be  accompanied  by  •: 
equivalent  charge  of  cations  to  maintain  electroneutrality .   It  is 
supposed  that,  in  the  absence  of  readily  available  cation  pools,  sol] 
neutralization  will  be  incomplete  and  that  dissolved  labile  monomeric 
aluminum  will  be  transported  from  the  soil  to  partially  counter  balance 
the  anionic  solutes  (Driscoll,  1985). 

The  ion  exchange  capacity  of  the  soil  appears  to  be  seasonally 
responsive  and  may  exert  greater  control  over  concentrations  of  alumina;:; 
than  solubility  reactions  would  predict  (Klock  and  LeFohn ,  1985). 
Differences  among  lakes  in  sensitivity  and  water  chemistry  are  largely 
the  result  of  variations  in  water  flow  through  soil  horizons  and    in  the 
i  ont  ict  of  the  water  with  surflcial  geology. 

There  are  several  mitigating  complexation  reactions  that  limit  the 
applicability  of  attributing  elevated  total  aluminum  concentrations  t  i 
acidification.   Aluminum  forms  complexes  with  sulfate,  fluoride,  organic 
carbon  compounds,  and  hydroxide  ions.   These  complexes  may  be  significant 
in  concentration  thereby  increasing  total  aluminum  solubility  independent 
of  the  rate  of  acid  loading.   Table  2  summarizes  the  aluminum  chemical 
equilibria  and  corresponding  thermodynamic  data  applied  in  this  study. 
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Table   2 


Thermodynamic  equilibrium  equations 
used  in  this  study  (DriscoJl,  1984) 


Equi libr ium 
Reaction  Constant 

Hydroxide  Ligands 

A13+-H20  =  A1(0H)2~-H~  1.03X10"5 

A13++2H20  =  Al(OH)2+-2H^  9.36xl0_11 

Al3++4Ho0  =  A1(0H.)"*4H*  6.93xl0"23 


Aluminum  Tr ihydroxide 

Al(0H)3(s)+3H'  =  A13+-3H?0 

Synthetic  Gibbsite  l.29x108 

Natural  Gibbsite  5.39x10° 

g 

Microcrystall ine  Gibbsite  2.24x10' 

Fluoride  Ligands 

Al3"-F~  =  A1F2~  1 .05x10' 

A13++2F~  =  A1F  "  5.77xl012 


A13+-3F   =  A1F„  ]  07X101 


A13'-4F"  =  A1F  !  5.37xl019 

4 

T+    -         ?-  ?0 

Al   +5F   =  A1F  c  8.33xl(T 

o 

A13+-6F"  =  A1F  3~  7.49X1020 

6 

Sulfate  Ligands 

A13  +  -S0  2~  =  A1S0/  1.63xl03 

4  4 

A13++2S0  2~  =  A1(S0  )  ~  1.29X105 

•-1  4  — 


3.0   METHODS  AND  PROCEDURES 

3.  1   Study  Area 

All  lakes  studied  are  in  Mount  Rainier  Nationa]  Park,  Washington, 
located  56  km  (35mi)  southeast  of  Tacoma.   Mount  Rainier  Park  covers  979 
Km2  (378  mi2)  surrounding  the  4,392  m  (14,410  ft)  peak.  Annua] 
precipitation  is  high,  151  to  254  cm  (60  to  100  in),  much  of  it  falling  as 
snow  at  this  high  elevation.   The  park  has  numerous  lakes  with  small  drain- 
age basins  characterized  by  shallow  soils  and  weathering  resistant  bed- 
rocks . 

Lakes  of  various  morphologies  were  sampled  from  all  sides  of  Mount 
Rainier,  with  elevations  ranging  from  1,200  to  over  2,000  m  (4000  to  6800 
ft)  as  shown  in  Figure  1.   Drainages  varied  in  size,  forest  type',  and 
forest  age.   Most  lakes  could  be  characterized  as  shallow,  (<  16  m)  with 
small  surface  area/watershed  ratios,  low  productivities  and  short  water 
residence  times.   These  lake  characteristics  are  indicative  of  sensitive 
lakes  (Driscoll,  1985;  Galloway  et  al . ,  1983;  and  Henricksen,  1980).   Based 
on  a  previous  study  (U.S.G.S.,  1986)  six  key  lakes  were  identified,  three 
of  which  had  high  total  aluminum  concentrations  (Golden,  Unnamed,  .<.  :<: 
Shriner),  and  were  matched  by  drainage  type  to  three  lakes  which  had  low 
total  aluminum  concentrations  (Mowich,  Crystal,  and  Allen).   Ten  other 
lakes  were  selected  to  provide  a  spectrum  of  drainage  and  lake  characteris- 
tics.  Table  3  summarizes  the  physical  characteristics  of  the  study  lakes. 


LAKE  EUNICE 


Q  MOWICH 
LAKE- 


GOLDEN 

LAKE 
o 


°LAKE   '.. 
ELEANOR 


UNNAMED 
LAKE 


,.  LAKE 
GEORGE 


°ALLEN  LAKE 


..  A*-'--"' 


REFLECTION-  "LOUISE 
LAKE       LAKE 


O        I        2       3      4      5 

I 1 L_l L_J     Mi. 


0     1     2    3   4     3 

1__1 1—1 I I    Km. 


FROZEN  LAKE 


SHADOW 

LAKE...  •  -ifCRYSTA^ 

>•?•-•'!■••'••  $■     \    LAKE 

P:   tipsoo* 

<~Z£'     LAKE:? 


FAN  LAKE 


oMARSH  LAKE 


Figure  1.   Study  Lakes  Mount  Rainier  National  Park 
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TARLF  3 


Physical  Characteristics  of  the  Lakes  Studied 


Relat  :\e- 
Max.   Sampling        Inflow    Outflow  Vegetation 
Lake      Altitude   Depth   Depth    Area   Observed   Observed   Density 
m        m       m       ha 


Golden 

1389 

6.  1 
r 

24.1 

7.3 

Y 

Y 

thick 

Unnamed 

1379 

3.7 

r 

11  .43 

0.8 

N 

Y 

thick 

Shriner 

1306 

6.1 
e 

2.8 

3.6 

I 

Y 

thin 

Crystal 

1777 

22.9 
r 

9.  i 

3.6 

I 

I 

modera te 

Allen 

1410 

7.6 
e 

7.6 

2.0 

N 

I 

modern te 

Mowich 

1289 

59 
m 

30.5 

49.6 

Y 

Y 

moder  i4>  e- 

Louise 

1400 

18.3 
e 

17.4 

6.9 

Y 

Y 

moderate 

Shadow 

1390 

4.6 
e 

4.9 

1  .2 

X 

Y 

thin 

George 

1290 

45.7 
r 

30.5 

1 3  .  3 

X 

X 

modera te 

Marsh 

1204 

6.  1 
e 

3.2 

1  .6 

X 

Y 

thick 

Tipsoo 

1620 

4.6 
r 

2.4 

2.6 

X 

Y 

thin 

Frozen 

2057 

5.5 
r 

6.1 

1.01 

X 

Y 

thin 

Eunice 

1632 

sh 

19.2 

4.9 

X 

Y 

moderate 

Reflection 

1482 

18.2 
V 

10.  1 

5.3 

X 

y 

thin 

Eleanor 

1512 

12.2 
e 

12.8 

7.9 

X 

y 

modera i c 

Fan 

1648 

sh 

0.5 

2.0 

X 

Y 

th  in 

r  =  report 

ed ,  unverified 

Y  =  yes 

e  =  estima 

ted 

X  =  no 

m  =  measured 
sh  =  shallow 
I  =  intermittent,  July  sampling  trip  only 
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3.2   Field  Sampling 
3.2.1   Data  Collection 

The  key  lakes  were  sampled  three  times  each  between  July  3  anil 
September  27,  1985.   The  ten  other  lakes  were  sampled  one  time  each  during 
the  same  period. 

Water  samples  were  collected  and  stored  in  1  liter  (1.0G  qt) 
polyethylene  bottles.   Prior  to  collection,  the  sample  bottles  were  washed 
with  1  N  nitric  acid,  thoroughly  rjnsed  with  deionized  water,  and  final- 
rinsed  on  site  with  water  from  the  lake  to  be  sampled.   Samples  were  stored 
on  ice  within  three  hours  of  collection  and  subsequently  refrigerated  at 
4  C  in  the  lab  trailer  and  home  lab  until  analysis. 

Water  samples  were  collected  at  two  depths  for  the  fifteen  spatially 
sampled  lakes  using  a  2  liter  (2.15  qt  )  Van  Dorn  bottle  water  sample!  .   Two 
1  liter  (1.06  qt)  "surface"  samples  were  collected  1  meter  (3.3  ft)  below 
the  surface.   A  1  liter  (1.06  qt)  "depth"  sample  was  collected  1  meter  (3.3 
ft)  above  the  lake  bottom  or  at  30.5  m  (100  ft),  whichever  was  less. 
Samples  were  obtained  from  the  region  of  maximum  depth  as  determined  by  a 
weighted  and  calibrated  drop  line.   In-situ  pH  was  determined  from  a 
subsample  of  the  collected  water  sample. 

Temperature  and  dissolved  oxygen  (DO)  profiles  were  determined  from  all 
spatially  sampled  lakes  at  0.76m  (2.5  ft)  or  1.5m  (5  ft)  intervals  to  a 
maximum  depth  of  15.24  meters  (50  ft)  using  a  polarographic  oxygen  meter 
(Y.S.I.  Model  54)  (Y.S.I.,  1980).   Prior   to  sampling,  the  instrument  was 
calibrated  on  site  for  altitude  (as  determined  from  U.S.G.S.  quadrangle 
maps)  and  for  ambient  temperature.   Equilibration  time  for  DO  and 
temperature  readings  was  three  to  five  minutes  or  greater. 

15 


All  samples  fit  Fan  Lake  and  the  August  sample  at.  Tipsoo 

lake 

were 

grab 

sampled  from  the  shore  at  a  de 

pth  of  0.3m.   Sampling  location 

was 

near 

the 

lake  outlet  and  was  selected  t 

o  avoid  areas  where  the  physica 

1  characti 

;r~ 

is  tics  differed  from  those  of 

the  general  littoral  zone.   Tern 

perat 

ure  and 

in-situ  pH  data  were  also  coll 

ected  at  the  sampling  location. 

All  in-situ  pH  measuremen 

ts  were  made  with  a  portable  di 

gital 

pH 

neter 

(Orion  Model  211)  and  a  combination  electrode  (Orion  Model  91 

-05). 

Th. 

>  pH 

meter  was  standardized  on  site 

at  pH  7  and  the  slope  was  cali 

brate 

d  at 

pH  4 

(Orion,  1981).   To  avoid  tempe 

rature  Induced  error,  the  pi!  of 

the 

buffi 

j 

and  samples  was  determined  at. 

their  respective  temperatures. 

To  1 

i  in  i  t 

the 

error  induced  by  C0o  exchange 

between  sample  water  anil  the  .it 

mosphere  ; 

3 

narrow  necked   125  ml  plastic 

bottle  was  used  for  in-situ  pH 

measurement. . 

Field  pH  was  determined  immedi 

ately  upon  retrieving  the  sampl 

e.  Th 

e  p!I 

and 

temperature  equilibration  time 

varied  from  five  to  fifteen  mi 

nut  e 

3.2.2   Lab  Trailer 

A  mobile  laboratory  trail 

er  provided  facilities  in  the  f 

i  e  1  d 

for 

equipment  preparation,  instrument  calibration,  and  initial  sample 

process  ing 

and  storage.   Aluminum  samples 

were  frud  i  ona  ted  into  I  abi  1  e 

and  a 

on  3  a! 

(organically  complexed)  components  and  fixed  for  later  analysis. 

3.3   Laboratory  Analyses 

All  water  samples  collect. 

ed  were  measured  for  laboratory 

equi 

libr; 

i  ted 

pH ,  conductivity,  alkalinity, 

major  metal  cations  (calcium,  magnesium, 

sodium,  potassium),  fluoride, 

total  aluminum,  chloride,  sulfa 

te,  t 

otal 

dissolved  organic  carbon  (DOC) 

and  dissolved  silica.   In  add! 

t  ion  , 

organically  complexed  aluminun 

was  determined  for  surface  sum 

10 

p  1 1". . 

3.3.1  Laboratory  pH 

Lab  equilibrated  pll  was  determined  in  conjunction  with  alkalin.it] 
titrations.   Three  digital  pH  meters  (Orion  Model  601-A)  with   combination 
electrodes  (Model  91-05)  were  used  (Orion.  1977A) .   To  ensure  consistency, 
results  were  verified  using  a  single  pH  meter  with  a  Ross  combination 
electrode  (Orion  Model  81-02).   The  pH  meters  were  calibrated  at.  pH  7,  and 
the  slope  was  adjusted  at  pH  4  at  ambient  temperature  immediately   before 
use  and  between  each  sample  measurement.   Samples  were  allowed  to  equili- 
brate at  ambient  temperature  prior  to  analysis. 

To  determine  pH .  30  ml  of  sample  water  was  placed  in  a  100  ml  plastic, 
beaker  and  stirred  magnetically.   An  insulating  pad  was  placed  between  the 
beaker  and  magnetic  stirrer  to  limit  heat  transfer.   Equilibration  of  pH  re 
quired  60-90  minutes,  and  the  data  were  recorded  after  the  pH  ren  M 
remained  stable  for  a  minimum  of  three  minutes. 

3.3.2  Conductivity 

Conductivity  as  umhos/cm  was  determined  using  a  cell-type  conductivity 
meter  (Lab  Line  Lectro-Mho-Meter ,  Model  MC-2  Mark  IV).   All  samples  were 
analyzed  for  conductivity  within  three  days  after  field  sampling.  Samp] 
were  either  allowed  to  equilibrate  at  standard  temperature  (25  C)  prior  to 
determination  or  the  data  were  corrected  to   25°C  using  linear  correlations 
developed  in  the  laboratory. 

3.3.3  Alkalinity 

Arid  neutralizing  capacity  ( [ANC]*- .carbonate) is  determined  by  the 
equivalents  of  acid  required  to  titrate  the  solution  to  a  known  equivalence 
po in  1  or  proton  reference  level  rind  is  the  parameter  measured  b\  the 
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titration.   In  carbonate  dominated  systems  alkalinity,  defined  as  [HC0o  ] 

2[C0^2~]  •-  [OH"]    [Hi  is  equal  to  [AXC]      Since  in  the  pH  range  3.5 

to  9.3,  [HCO  ]  is  much  greater  in  concentration  than  the  othe 

constituents  its  molar  concentration  can  be  assumed  to  equal  alkalinity. 

There  was  no  significant  difference  in  analytical  results  by  choosing  one 

parameter  over  the  other. 

Alkalinity  ([ANC]_    for  the  carbonate  system)  was  determined  by 

potentiometrically  titrating  a  50  ml  water  sample  with  0.02N  HC1  using  an 

adjustable  micropipet  tor .   Samples  were  stirred  at  low  speed  to  limit  C0o 

stripping  during  the  titration.   A  minimum  of  five  and  up  to  fifteen  titi 

Lion  data  points  (total  equivalents  of  acid  vs.  pH)  were  regressed  to 

determine  the  [ANC]-__  equivalence  point  using  linear  Gran  functions. 

(Gran,  1950;  Stumm  and  Morgan,  1931).   All  linearly  regressed  tltratii 

data  points  were  within  the  pH  range  3 .  5  Lo  3  to  limit,  both  weak  basi 

hydrogen  activity  errors  inherent  in  this  method  (Driscoll  and  Bisogni, 

1982).   Analyses  were  normally  performed  in  duplicate,  but  La  triplicate  if 

the  relative  error  between  duplicates  exceeded  10°,..  Results  typically 

agreed  within  5%   of  the  mean.  Standard  solutions  of  40  meq/I.  [ANCl-  _  were 

i  =u 

titrated  to  verify  acid  normality.   [AXC"     deterirunal  :<o  '■■     e  i  up]  >t< 
within  five  days  after  field  sampling. 

3.3.4   Major  Cation  Metals 

2+     2      - 
Analyses  for  major  cation  metals  (Ca   ,  Mg   ,  Na  ,  K  )  were  performed 

using  an  atomic  absorption  spectrophotometer  (Perkin-F.lmer  Model  360).  Th< 

ten  second  integrating  cycle  was  used  for  quantification.   Lanthanum  (2  "" 

in  50°..  PfCl)  was  added  to  the  samples  .xml    standards  to  eliminate 

interferences  asso<  Lated  with  the  ionization  of  alkaline  ea  'I.-.,    di  solved 
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silica,  and  aluminum.   A  double  glass  distilled  water  blank  was  carried 
through  all  analyses  to  provide  a  basis  for  background  correction.   Samples 
were  pretreated  by  filtration  through  a  C.l  urn  membrane  filter  and  there- 
fore reflect  dissolved  concentrations.   Detection  limits  of  130  ug/L 
calcium,  15  ug/L  sodium,  13.5  ug/L  potassium,  and  6  ug/L  magnesium  were 
obtained  by  this  method. 

3.3.5   Fluoride 

Fluoride;  analyses  were  performed  using  a  fluoride  specific  ion 
electrode  (Orion  Model  96-00)  and  an  Orion  Ionanalyzer  (Model  801)  capab] 
of  reading  to  0.1  millivolts.   To  eliminate  error  induced  by  changes  in 
fluoride  activity,  50  ml  of  low  level  total  ionic  strength  adjusting  buffe! 
(TISAB)  was  combined  with  "0  ml  of  sample  (Orion,  1077B) .   Fluoride 
response  was  linear  with  concentration  down  to  100  ug/L  one  curvilinear 
from  100  ug/L  to  the  detection  limit  of  1  ug/L.   At  the  low  concentrations 
encountered  in  most  lake  samples,  electrode  equilibration   required  betweei 
20  and  45  minutes. 


3.3.6.1.   Aluminum  Analyses 

Aluminum  was  analyzed  by  an  adaptation  of  Barnes  (1075)  method  (La 
Zerte,  1984).   Aluminum  in  a  200  ml  filtered  sample  was  complexed  with  2  ml 
of  5°o  8-hydroxyquinol  ine  and  buffered  at  pH  8.3.   The  resulting  complex  was 
then  extracted  into  10  ml  of  methyl isobutylke tone  (MIBK)  by  mixing  foi  10 
seconds  followed  by  3  to  5  minutes  separal  Ion   The  lower  aqueous  laye   was 
then  drawn  off  using  a  separatory  funnel  and  discarded.   The  organic  laye; 


19 


was  stored  at  4°C  in   dark  brown  glass  vials  until  analysis.   Aluminum 
determinations  were  made   using  a  spectrophotometer  (Heath  Series  700)  at 
003.01  nm  using  1  cm  quartz  colls. 

Reaction  time  of  aluminum  extraction  was  minimized  by  the  ten  secoiK'. 
mixing  time.  Therefore,  only  the  monomeric  (free  or  complexed)  aluminum 
species  that  can  be  considered  as  equilibrium  species  were  removed  (Barnes, 
1975;   Smith  and  Hem,  1972). 

Surface  samples  were  fractionated  between  the  operationally  defined 
labile  (inorganic)  and  nonlabile  (organically-complexed)  aluminum  by  passing 
an  aliquot  throng!)  a  column  of  strongly  acidic  cation  exchange  resin 
(Amberlite  C.G.-120).   The  column  was  regenerated  with  10  "  y.   N'sCl    '• 
each  sample  separation  to  maintain  the  resin  in  the  sodium  form.   The  flow 
rate  was  maintained  above  25  ml  'min  to  limit  the  disruption  of  labile 
complexes  within  the  resin  (Driscoll,  1004).   Organic  ligands  form  stri 
complexes  with  aluminum  and  are  presumed  to  be  anionically  or  nonionicalli 
charged,  thus  effectively  competing  with  the  exchange  resin.   The  free  and 
inorganically  complexed  aluminum  are  readily  removed  from  solution  by  the 
polar  resin.   The  pH  change  associated  with  the  cation  exchange  treatment  is 
sufficient  to  disrupt  the  cationically  charged  inorganii  ligands  (A1(0H), 
but  will  only  affect  the  nonlabile  monomeric  aluminum  to  a  limited  extent 
(Driscoll,  1984).   The  nonlabile  aluminum  is  measured  directly  in  the  column 
eluate.   The  difference  between  total  monomeric  aluminum  and  nonlabile 
aluminum  is  termed  inorganic  aluminum.   This  fraction  includes  free  aluminum 
as  well  as  hydroxide,  fluoride,  and  sulfate  complexes  of   aluminum. 

-.0.0.2   Storage  Procedure 

Aluminum  storage  procedures  were  tested  for-  accuracy  and  t.l;::e  indiu  ed 
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error.   The  standard  curve  was  found  to  be  linear  down  to  7  ug  L  and 
logarithmic  below  7  ug 'L  to  a  detection  limit  of  0.35  ug/L .   Eighty-one  per 
cent  of  the  samples  fell  on  the  nonlinear  part  of  the  curve.   The  standard 
error  of  0.28  ug/L  on  the  lower  portion  of  the  curve  resulted  in  a  9"°,, 
confidence  interval  of  less  than  20°o  of  the  mean. 

Field-fixed  samples  that  were  stored  three  months,  then  reanalyzed, 
averaged  ll°i  higher  total  aluminum  concentration  than  the  original 
determination.   This  change  is  not  considered  important  when  compared  to  the 
original  confidence  interval  for  aluminum  determination  of  20°..  of  the  mean. 
Aluminum  fixed  three  months  afte]  sample  collection      iged  26"  !owu:  thar. 
the  original  sample  (fixed  three  to  four  hours  after  collection)  and  had 
increase  in  the   confidence  interval  to  36°u  of  the  mean.   Apparently  some 
monomeric  aluminum  was  "lost".   Aluminum  samples  for  the  July  trip  were 
fixed  upon  returning  to  tin;  laboratory  one  to  three  weeks  afte   •  ■  ...  1 .!  , 
may  therefore  be  slightly  underestimated.   Recent  research  indicates  the 
loss  of  aluminum  is  associated  with  the  organically  complexed  fraction 
(Sullivan  et.  al..  1986).   Equilibrium  may  also  have  shifted  in  favor  of  the 
polynuclear  aluminum  species  not  determined  by  this  method  (Barnes.  1973 
Sullivan  et.  al . .  1986) . 

3.3.7  Chloride 

Chloride  ion  concentration  was  determined  potent jometrical ly  in  accord 
ance  with  Standard  Methods  No.  407  B  (A.P.H.A..  1983). 

3.3.8  Sulfate 


Sulfate  concentration  was  determined  using  a  variation  of  the 
turbidimetric  method  (A.P.H.A.,  1933)  as  described  by  N'elson  and  Delwii 
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(1983).   A  turbidity  meter  (Hach  Mode]  2100-A)  was  used  for  turbidity 
measurement . 

The  turbidimetric  method  involves  the  iddition  of  barium  chloride  to 
the  sample,  causing  the  precipitation  of  the  less  soluble  barium  sulfate 
The  conditioning  agent  described  in  Standard  Methods,  normally  used  to  as 
sure  crystals  of  uniform  size,  was  not  used  in  this  procedure  because  of  its 
associated  background  turbidity.   Eliminating  the  background  turbidity 
increased  the  precision,  lowered  the  detection  limit,  and  increased  the  re- 
producibility of  this  method  (Nelson  and  Delwiche,  1983).   The  standard 
curve  however,  is  exponential  rather  than  linear  due  to  the  large   be  [ u n 
sulfate  crystals  at  higher  concentrations.   As  a  result  .  numerous  .    lards 

were  required  to  define  the  curve  in  the  sampling  range.   The  detection 

2~ 

limit  was  <0.05  mg/L  as  SO  "  . 

4 

3.3.9  Total  Dissolved  Organic  Carbon  (DOC) 

Samples  were  filtered  through  pre-combus ted  glass  fiber  filters 
(Whatman  934-AH).   DOC  was  determined  instrumental ly  by  CO   infrared 
absorbance  of  wet  oxidized  dissolved  organic  carbon  (Oceanography 
International  Model  0524AA  Carbon  Analyzer)  and  prepared  in  accordance  ;<  :  '  ' 
the  Oceanography  International  Carbon  Analyzer  Instrument  Manual.   All 
ampules  used  for  samples  or  standards  were  precombusted  at  300  C  for  12 
hours.   Samples  were  analyzed  for  DOC  with  storage  time  ranging  from  three 
to  five  months.   The  detection  limit  was  0.5  mg/L  as  carbon  with  a  standard 
error  equal  to  17°o  of  the  mean. 


1.3.  .10   Silica 

Dissolved  silica  concentrations  were  determined  spei  trophol  ■■■:■'■ 
Bausch  and  Lomb  Spectronic  3,°,)  <u  315  nm  using  the  he  tempo  ly  bi  ie  method 


(Standard  Methods  No.  425  D.  A.P.H.A.,  1983).   The  detection  limit  was  0.02 
mg/L  as  Si0o. 

i . 0   RESULTS 

The  lakes  sampled  in  this  scudy  may  be  characterized  as  weakly 
buffered,  with  an  extremely  dilute  ionic  composition  (conductivity  range 
3.6  to  56.3  uhmos/cm) .   Calcium  was  the  dominant  cation,  constituting  76.4% 
of  the  total  cations  on  a  equivalent  basis,  followed  by  sodium,  16.26%; 
magnesium,  6.02%;  and  potassium,  1.32V   Bicarbonate  dominated  the  anions, 
accounting  for  80.01%  on  an  equivalent  basis  followed  by  sulfate,  3.3,12%; 
and  chloride,  6.47%.   Composited  lake  surface  water  chemistry  data  ire  sum- 
marized in  Table  4,  with  individual  lake  water  chemistry  data  summarized  in 
Appendix  3. 

J . 1   Limnology 

Limnological  data  is  presented  in  Appendix  1.   Lakes  appeared  to  reach 
a  thermal  maximum  in  July  or  August  and  decrease  in  temperature  by 
September.   Stratification  is  more  pronounced  in  July  than  September.   Fall 
turnover  appears  to  have  been  completed  by  September,  with  only  Golden  Lake 
remaining  stratified  at  the  end  of  the  summer.   Dissolved  oxygen 
concentrations  generally  increased  as  lake  temperatures  decreased. 

Temporally  sampled  lakes  had  lowest  alkalinities  in  July,  increased  to 
a  maximum  in  August  and  generally  showed  a  decrease  by  late  September  as 
shown  in  Figure  2.   The  most  significant  variation  in  alkalinity  occurred 
in  small  volume  to  watershed  area  lakes.   The  higher  volume  to  watershed 
area  ratios  of  Golden  and  Mowich  lakes  result  in  longer  hydraulic  residence 
times  which  dampen  seasonal  events  (Welch  et  al.,  1985). 
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Composited  Water  Chemistry  Summary 
for  Lake  Surface  Samples 


Parameter 

Units 

Max . 

M  i  n . 

Ave, 

Alkalinity 

ueq/L 

404 

17 

106 

[so/"] 

ueq/L 

99 

0 

1  ] 

[Ca2~] 

ueq  I. 

403 

10 

93 

[Mg2+] 

ueq/L 

27 

0.7 

6.4 

[NaT] 

ueq  I. 

32 

:o 

19 

[K+] 

ueq/L 

0.5 

I  .  4 

Anion/cation 

ratio 

1  .26 

0.S3 

0.93 

Conductivity 

umbo /cm 

36.3 

3.6 

16.2 

[HCO  " ]  '[cation; 

5] 

1  .0 

0.76 

0  36 

pH  (field) 

7.3 

5.3 

6 .  3  ~ 

pH  (lab) 

7.5 

3.6 

6.6 

DOC 

mg.'L 

3.09 

0.5 

-«   O 

* 

Lakes  sampled  multiple  times  were  averaged  prior  to  computing 

the  overall  summary.   All  values  are  sea  salt  corrected. 


The  drop  in  alkalinity  seen  in  September  for  most  temporally  sampled  lakes 
is  often  associated  with  fall  turnover  (Lazerte,  1984). 

There  was  a  definite  increase  in  alkalinity  and  a  generally  greater  ioi 
concentration  with  depth  in  all  lakes.   The  deeper  samples  (generally  below 
the  thermocline)  exhibited  bettei  correlation  coefficients  between  ion 
concentrations  and  conductivity  than  did  their  surface  sample  counterparts . 
indicating  a  greater  dependence  of  ionic  composition  on  mineral  dissolution 
These  profiles  are  generally  consistent  with  the  physical  behavior  of 
dimictic  lakes  and  suggest  that  during  the  summer,  water  more  typical  of 
groundwater,  without  precipitation  dilution,  accumulates  in  the  hypo].! hid 
(Welch  et  al . ,  1985) . 

4.2  pH 

The  range  of  pH  values  determined  for  lakes  in  this  study  was  3  3 
7.4  in  the  field  and  3.6  -  7.6  for  lab  equilibrated  samples.   Lab- 
equilibrated  pH  averaged  0.08  pH  units  higher  than  field  pH  with  a  standard 
deviation  of  0.012  pH  units.   Empirically  derived  equations  correlating  the 
parameters  measured  in  this  study  to  lab-equilibrated  pH  are  listed  in 
equations  1-16  (App.  2).   Parameters  are  sea  salt  corrected  and  expressed 
in  ueq/L  except  where  noted. 

4.3  Conductivity 

The  average  conductivity,  for  all  samples,  of  22.7  umhos/cm  reflects 
this  dilute  nature  of  the  lakes  studied.   Conductivity  ranged  between  3.4 

and  83.0  umhos/cm  with  a  standard  deviation  of  19.8  umhos/cm.   Conductivity 

2 

is  strongly  correlated  with  ionic  strength  (r    0.93)  and  therefore 

provides  an  indirect  measurement  of  the  water's  overall  ionic  composition. 
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A  strong  correlation  between  conductivity  and  individual  ionic  conc- 
entrations indicates  weathering  reaction  genesis.   OnJy  tot . j  1  aluminum  and 

2 

chloride  showed  no  significant  correlation  with  conductivity  (r   ;  0.30). 

Correlation  coefficients  for  all  of  the  ion  concentration  to  conductivity 
equations  increased  with  lake  depth. 

The  correlation  between  ionic  strength  (mole/L)  and  conductivity 
(umhos/cm)  for  this  study  (T  =  -3.09  +  1. 44x10  v)  x  cond;  eq  .  23)  is  similar 
to  both  the  equation  determined  for  the  Oregon  Cascade  lakes  (I  =  1.28x10 
x  cond;  Nelson  and  Delwiche,  1983)  and  the  equation  cited  by  Snoeyink  and 
Jenkins  (1900)  for  surface  waters  in  general  (I    1.6x10  "  x  cond.). 
Activity  coefficients  of  0.93  and  C.933  fur  monovalent  and  divalenl  ions, 
respectively,  were  calculated  by  the  Debye-Huckel  limiting  law  for  the  mean 
lake  ionic  strength.   Errors  incurred  by  using  concentrations  rather  than 
activities  are  small  and  generally  within  the  analytical  uncertainty  of  Ilia 
data . 

1.4   Alkalinity 


The  mean  alkalinity  for  all  samples  was  162  ueq/L  with  a  range  of  10.3 
to  693.5  ueq/L  and  a  standard  deviation  of  163  ueq'L.   For  lake  surface 
samples  tin;  median  and  average  alkalinjties  ( [HCO   j)  were  62.5  ueq/L  and 
106,  ueq/L  respectively.   These  results  are  similar  to  those  found  by  Welch 
and  Chamberlain  (1981)  for  18  lakes  north  of  Mount  Rainier  in  Western 
Washington  state  (median  =  73  ueq/L  and  average  113  ueq/L). 

Alkalinity  correlations  with  measured  parameters  are  listed  in 
equations  29-41  (App.  2).   Alkalinity  as  a  function  of  lab  pH  and  of  field 
pH  are  illustrated  in  Figures  3  and  4,  respectively.   Lake  data  is  plotted 
and  identified  by  sample  depth.   The  theoretical  alkalinity  vs.  pi!  equation 


(I/ben)  A4iUi|D>i|D 

28 


(|/ben)  A4iui|D>||D 
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for  mean  sea  level  prn   =  jO  '    atm.  is  mciuded  fur  comparison.   The 

2 
data  more  closely  correlate  with  the  theoretical  equation  (eq.  16,  Table  1} 

when  lab-equilibrated  pH  is  used.   Points  falling  above  the  theoretical 
line  indicate  either  equilibration  at  a  higher  C0„  partial  pressure,  such 
as  from  groundwater,  or  the  presence  of  an  additional  weak  buffering  sys- 
tem.  These  relationships  are  discussed  in  the  section  entitled  "Buffering 
System"  (Sec .  5.3) . 

4.5  Major  Cation  Metals 
4.3.1   Calcium 

Calcium  concentrations  averaged  2.3  mg/L,  ranged  between  0.2  and 

3.6  mg/L  and  had  a  standard  deviation  of  2.9  mg/L.   Calcium  was  the 
dominant  cation  in  all  but  a  few  exceptionally  dilute  lakes  where  sodium 
dominated.   Strong  correlations  were  found  between  calcium  and  alkalinity, 
conductivity,  and  lab  equilibrated  pH  (r~  >  0.81).   Calcium  increased  with 
depth  for  almost  all  lakes  and  surface  concentrations  increased  over  t.he 
summer  sampling  period.   Of  the  cations,  calcium  most  closely  reflected  the 
temporal  variations  of  alkalinity. 

A  comparison  of  alkalinity  to  calcium  ratios  for  several  geographi 
areas  is  presented  in  Table  3.   This  relationship  is  remarkably  consistent 
for  Washington  Cascade  lakes.   Logan  et  al.  (1933)  concluded  that,  assuming 
similar  geology,  the  similarity  in  equations  indicates  either  no 
acidification  or  consistent  magnitudes  of  acidification  among  Washingtion 
Cascade  Lakos . 

Similar  Lo  the  more  rigorous  eiectroneutralit.y  filiation,  the  intercept 
term  of  these  simplified  equations  may  be  used  with  caution  to  indicate 
strong  acid  concentrations.   The  greater  the  effecl  of  strong  acid,  !':••  no 
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negative  the  intercept  term.  Apparently  there  is  a  minima]  effect  of  strong 
acid  deposition  on  Washington  and  Oregon  lakes.   The  significantly  lower 
intercept  term  for  the  described  regression  applied  to  Adirondack  lakes 
(New  York)  data  indicates  the  greater  effect  of  acids  and  probable  anthropo- 
genic influences  on  their  water  chemistry  (cited  in  Nelson  and  Delwiche, 
1983) . 

Table  5 
Alkalinity  vs.  Calcium  for  Several  Geographical  Are. 


Region 


Equat i  on 

(y=i0-b(x)) 


Correlation 

Coeff  .  ,  r 


Mt .  Rainier 


[A Ik]  =  1.83  -  1.09[Ca 


9 


Central  Washington    [Alk]=-5 . 77-1 . 07[Ca   ] 
Cascades 


North  Washington 
Cascades 

Oregon  Cascades 
1983 


[Alk]=2.R7-0.995[Ca-_] 


[Alk]=-3.70+2.40[Ca   ] 


Experimental  Lakes    [Alk]=-32 .0-1 . 42[Ca   ] 
Western  Ontario 


Northern  Norway 


2+. 


[Alk]=-29.0+1.32[Ca~    ] 


?.  +  . 


Adirondack  lakes      [Alk] =-94 . 9+0 . 94[Ca~  ] 


0  .  97 
0  90 


0.88 


California  Cascades   [Alk]=12. 775+0. 8518[Ca]    0.97 
1983 


0.7G 

0  .  83 

N/A 


Source 


his  3 1 u  i 


•  .  ■   • 


1      1  o 


Logan  et.  a.l  .  ,  1032 

Nelson  and  Delwiche 

Me  lack  et:  al.  , 

Henr i  cksen .  1979 

Henri  cksen ,  1 979 
Driscoll,  1980 


Data  not  sea  salt  corrected. 
N/A  =  not  available. 


4.5.2  Magnes ium 

Magnesium  concentrations  for  ali  samples  averaged  0.31  mg/L,  ranged 
between  0.02  and  1.17  mg/L,  and  had  a  standard  deviation  of  0.30  mg/L. 

o 

Magnesium  was  found  to  be  closely  correlated  with  alkalinity  (r   -  0.92), 

2  ? 

and  conductivity,  (r   =  0.93)  and  have  a  fair  correlation  with  pH  (r 

0.73).   Correlation  coefficients  were  less  than  those  for  calcium  but 

greater  than  those  for  sodium  and  potassium. 

4.5.3  Sodium 

Sodium  concentrations  for  all  samples  averageu1  0.70  ug,  L,  had    .  mg< 
of  0.24  to  2.24  mg/L  and  a  standard  deviation  of  0.39  mg/L.   Sodium  showed  a 
weak  correlation  with  pH  (r   =  0.31)  and  fair  correlations  with  both 
conductivity  (r  =  0.67)  and  alkalinity  ( r2  -  0.73).   Stepwise  nultiple 
regression  showed  a  significant  correlation  increase  when  sodium  was 
included  as  a  component  relative  to  alkalinity.   Therefore  a  significant 
amount  of  sodium  appears  to  be  derived  from  mineral  dissolution  with  bicar- 
bonate as  a  counterion. 

4.5.4  Potassium 

Potassium  concentration  for  all  samples  averaged  0.10  mg/L,  had  a  range 

between  0.02  to  0.71  mg/L  and  a  standard  deviation  of  0.14  mg/L.   Potassium 

2 
showed  only  weak  correlations  with  conductivity  (r   =  0.43)  and  alkalinity 

2  ? 

(r   =  0.53)  and  a  low  correlation  with  pH  (r   =  0.22).   Similar  results  were 

reported  for  the  Oregon  Cascades  (Nelson  and  Delwiche,  1983).   Despite  its 

relatively  low  concentration  ( 1  . 3°n  of  cations)  stepwise  multiple  regression 

shows  potassium  to  be  the  most,  significant  counterion  for  alkalinity,  afte; 

calcium.   This  indicates  Lhal  potassium  is  the  product  of  weathering 
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reactions  . 

The  poor  correlation  coefficients  for  potassium  may  also  reflect  error 
induced  by  incorrect  sea  salt  estimates.   Potassium  in  aerosols  is  greatly 
influenced  by  terrestrial  dusts  which  are  not  accounted  for  by  sea  spray 
corrections,  as  applied.   Due  to  the  low  potassium  concentrations,  this 
error  may  be  significant. 

4.5.5   Total  Cations 

Total  cations  (not  sea  salt,  corrected,  surface  and  depth  samples) 
averaged  200.4  ueq/L,  ranged  between  21.8  and  741.8  ueq/L  and  had  a  standard 
deviation  of  185  ueq/L.   On  a  sea  salt  corrected  equivalent  basis  the 
Ca:Na:Mg:K  ratio  is  80.0:13.2:5.5:1.3.   Multiple  regression  showed  that  all 
major  cations  are  positively  correlated  with  bicarbonate  .;::  !  increase  the 
correlation  coefficient  in  the  descending  order  of  significance.  Ca>K  N'a  Mg. 
This  indicates  that  the  dominant  cations  are  all  counterions  for  bicarbonate 
in  weathering  reactions  (Zimmerman,  1982). 

High  correlations  were  found  between  total  cations  and  alkalinity  (r   = 
0.99),  conductivity  (r2  -  0.99),  and  pH  (r2  =  0.77). 

In  acidified  Adirondack  lakes  hydrogen  and  aluminum  account  for 
average  of  23°o  of  the  total  cation  concentration  (Driscoll,  1980).   The 
relative  cation  contribution  of  hydrogen  and  aluminum  ions  was  negligible  in 
this  study. 

4.6   Sulfate 

Sulfate  averaged  1.09  mg/L,  with  a  range  of  0.1  to  6.2  mg/L  and  a 
standard  deviation  of  1.5  mg/L.   Sulfate  bar!  c;  weak  positive  correlatio 

o 

with  conductivity  (r    0.46).  and  low  correlations  with  pH  ,\wi\   alkalinity. 
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These  relationships  are  greatly  affected  by  apparent  indigenous  sulfate 

sources  at  Allen  Lake.   When  Allen  Lake  data  is  disregarded,  sulfate 

2 
correlation  with  conductivity  is  poor  (r  <0.15), 

Sulfate  concentrations  did  not  increase  with  lake  depth.   The  lack  of 

conductivity,  alkalinity,  or  spatial  correlations  for  sulfate  indicates  that 

weathering  reactions  are  not  a  significant  source  of  sulfate.  Local  sources 

of  sulfate  resulted  in  a  negative  correlation  between  pH  and  [SO  ~  ]  in  the 

•4 

Alpine  Lakes  Wilderness  Area  of  Washington.   Similarly,  when  local  sources 
were  disregarded,  no  correlation  of  sulfate;  with  pH  was  evident  (Welch    ) 
Chamberlain,  1981).   No  significant  correlation  was  found  between  [SO,"  ] 
and  other  parameters  in  Oregon  Cascade  lakes  (Nelson  and  Delwiche,  198.°.) 

4.7  Chloride 

Chloride  averaged  0.42  mg/L  with  a  range  of  0.0"  to  1.12  ng  T  ,   ' 
standard  deviation  of  0.22  mg/L.   Chloride  accounted  for  an  average  of  6.5°. 
of  the  total  anions  on  an  equivalent  basis  and  showed  little  variation  in 
concentration.   No  significant  correlation  existed  between  chloride  and  p!-I . 
conductivity,  or  alkalinity  (r   <  0.22).  Sea  spray  appears  to  be  the 
dominant  source  of  chloride  in  the  study  lakes.   Sea  salts  (based  on  [CI  ]) 
accounted  for  11.535  of  total  cations  and  ll.l°o  of  the  anions  (equivalent 
basis)  in  the  study  lakes. 

4.8  Fluoride 


Fluoride  concentration  averaged  15.2  ug/L  with  a  range  of  6.5  to  40 
ug/L  and  a  standard  deviation  of  3.5  ug/L.   Fluoride  concentration  was 
positively  correlated  with  condud  ivity  >*n<}   alkalinity,  but  not  pH. 
Driscoll  (1984)  suggested  that  aluminum  complexation  controls  fluoi  Lde 


concentrations  in  Adirondack  lakes.   At  the  low  pH  of  Frozen  Lake,  nearly 
30°0  of  the  fluoride  present  was  estimated  to  be  complexed  with  aluminum. 

However,  at  the  pH  ranges  generally  encountered  in  this  study,  alumim 

fluoride  complexes  averaged  less  than  0.02°6  of  total  fluoride.  Average? 
fluoride  concentrations  of  15  ug/L  are  similar  to  the  11.8  ug/L  found  in 
Oregon  Cascade  lakes  (Nelson  and  Delwj'che,  1983)  and  are  probably  background 
levels  . 

4.9   Aluminum 


Aluminum  concentration  for  all  samples  averaged  14.8  ug/L,  ranged 
between  0.7  ug/L  and  90  ug/L  with  a  standard  deviation  of  25  ug/L.   For 
surface  samples  that  were  fractionated,  total  aluminum  averaged  13.2  ug/L, 
nonlabile  (organically  complexed)  aluminum  averaged  2.7  ug  L,  and  labile 

aluminum  averaged  10.6  ug/L.   Total  aluminum  was  only  weakly  correl  ited 

2 
with  alkalinity,  conductivity,  or  pH  (r  <  0.35).   Similar  or  improved 

correlations  with  aluminum  were  found  for  silica,  DOC,  and  field  p!i. 

Multiple  regression  indicated  that  over  65°i  of  the  variation  in  aluminum 

concentration  could  be  accounted  for  by  variation  in  DOC  and  field  pH. 

This  relationship  remained  consistent;  fin-  surface  and  depth  samples  and  was 

not  improved  by  the  addition  of  any  other  parameters. 

Surface  aluminum  samples  were  fractionated  into  organic  and  labile 

aluminum  components.   The  organically  complexed  aluminum  showed  a  good 

2  2 

correlation  with  DOC  (r   =  0.79)  and  was  independent  (r   <  0.25)  of  any 

other  parameter.   The  labile  (inorganic)  aluminum  component  showed  a  weak 

2  ° 

correlation  with  DOC  (r~  =  0.44)  ami  a  still  weaker  correlat  ion  with  pli  (r 

-  0.37). 


4.9.1   Aluminum  Speciation 

The  water  chemistry  data  for  aluminum  were  examined  with  respect  to 
theoretical  thermodynamic  and  mass  balance  relationships  as  previously 
summarized  (Table  2)  using  the  MICROQL  computer  program  developed  by 
Westall  (1979).   In  general,  fluoride  and  sulfate  complexes  with  aluminum 
were  negligible.   Fluoride  complexes  often  become  the  dominant  form  of 
aluminum  below  pH  5.5  (Driscoll  et  al.,  1984).   Only  Frozen  Lake  and  Colder: 
Lake's  July  sample  were  predicted  to  have  significant  aluminum  fluoride 
complexat ion .   The  predominant  aluminum  liganu  was  hydroxide  under  the 

conditions  encountered  in  this  study. 

3^ 
Free  aluminum  concentrations  ( [Al   ])  calculated  for  the  observed  data 

when  compared  to  theoretical  solubilities  for  aluminum  solid  phases  appear 

to  be  controlled  by  an  aluminum  trihydroxide  phase  similar  to  natural 

gibbsite  (log  k'sQ  =  -8.88)  and  less  soluble  than  microcrysta.il  .hie  gjbhslft 

(log   K    =  -9.35).   Calculated  free  iunic  aluminum  is  plotted  versus  pH 

and  showing  theoretical  solubility  controls  in  Figure  5.   The  solubility 

constant  determined  by  regression  for  the  lakes  in  this  study  (log  K   =  - 

8.86)  is  similar  to  the  constants,  log  K   =  -9.1  and  -8.49   empirically 

SO  f  -j 

determined  for  lake  systems  by  Cosby  et  al..  (1035)  and  Driscoll  as  ' 
Bisogni  (1982),  respectively. 

Figure  6  superimposes  the  labile  (total  inorganic)  aluminum  cohcentra 
tions  on  an  aluminum  solubility  diagram  in  equilibrium  with  macrocrystalline 

gibbsite.   Only  one  point  (Golden  Lake,  July)  would  theoretically  indicate 

3- 
elevated  levels  of  aluminum  ion  [AT  ]  due  to  pH  values  below  5.5.   Elevated 

levels  of  labile  aluminum  above  pi!  6  (Shriner  and  Tipsoo  Lakes)  are 

attributable  to  Al(OH)    complex  formation  as  pH  increases. 
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4.9.2  Temporal  Variation  of  Aluminum 

Aluminum  concentrations  were  lowest  in  August.  The  increases  in  both 

early  summer  and  fall  appear  to  be  largely  caused  by  increases  of 

organically  complexed  aluminum  resulting  from  increased  levels  of  DOC 

(Figure  7).   When  determined  seasonally,  the  apparent  solubility  constant 

decreases  from  log  K   =  -8.80  to  -9.04  from  early  summer  to  fall. 

so  •* 

Although  the  difference  between  these  two  values  is  probably  not 
significant,  it  does  indicate  a  trend  of  decreasing  aluminum  solubility 
over  the  summer-.   Driscoll  (1985)  attributed  similar  results  to  sp;  lug 
snowmelt  events  releasing  a  peak  flow  of  water  unlikely  to  be  in 
equilibrium  with  any  controlling  aluminum  solid  phase. 


10   Total  Dissolved  Organic  Carbon  (DOC) 

The  average  total  DOC  for  all  samples  was  2.2  mg/L  as  carlo;  with 
range  of  0.005  mg/L  to  24.8  mg/L  and  a  standard  deviation  of  4.9  mg/L.   The 
average  lake  surface  total  DOC  concentration  of  1.31  mg/L  as  carbon 
indicates  the  unproductive  nature  of  the  study  lakes.   A  similar  range  (0.1- 
3.6  mg/L  DOC)  was  found  for  six  Oregon  Cascade  lakes;  however,  four  of  these 
lakes  had  concentrations  below  th<j  detection  limit  of  0.1  mg  I.  (Nelson  and 
Delwiche,  1983).   The  DOC  concentrations  were  similar  to  or  below  levels 
reported  for  other  dilute  systems  (Elzerman,  1983).   For  surface  waters, 
levels  of  DOC  are  independent  of  field  pH  or  any  parameter  other  than 
organic  aluminum. 

Levels  of  DOC  were  an  order  of  magnitude  higher  in  the  Golden  :.ak<j 
depth  sample  than  those  encountered  in  the  other  samples  depicted  In  Figi  n 

7,  indical  ing  either  high  lake  produi  '  '  v  :  t  y  or  i haps  local  producl  ion  In 

tluj  soil  ar.d  runoff  into  the  lake. 
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4.11   Silica 


The  average  silica  concentration  for  all  samples  was  4.5  ma/l.    is  SID, 


c       •"»   /? 


Silica  exhibited  moderate  correlations  with  conductivity  and  aluminum  (r~  < 
0.60)  and  a  weak  correlation  with  pH  (r   =  0.25).   Silica  either  remained 
stable  or  showed  a  general  decrease  in  concentration  through  July  to  Septem- 
ber, a  pattern  attributed  in  other  lake  systems  to  uptake  by  diatom: 
(Driscoll,  1930;  Wetzel,  1975).   The  low  silica  concentrations  and  low 
correlations  between  silica  and  conductivity  or  alkalinity  obtained  fri 
surface  samples  indicate  either  unreactive  bedrock  or  that  much  of  the 
water  enters  predominantly  as  precipitation.   Such  conclusions  should  be 
made  with  caution  due  to  the  biological  mediation  of  the  silica  cycle 

(Driscoll,  1980).   Silica  correlations  with  conductivity  (  ~    n  05)  and 

p 

alkalinity  (r  =  0.07)  improved  with  lake  depth,  ind.icat  ing  sil  )c   i's 

vation  from  mineral  dissolution. 

4.  12   Buffer  System 

Buffer  capacity  may  be  defined  as  the  number  of  equivalents'  liter  of 
strong  acid  which  when  added  to  a  solution  causes  a  unit  c)  in  ;.:" 

Figure  8  presents  titration  curves  for  five  lakes  with  varying  alkalinJUes 
sampled  in  this  study.   An  experimentally  determined  buffer  intensity 
diagram  corresponds  directly  to  the  slope  of  these  titration  curves  and  was 
determined  from  discrete  increments  of  acid  added  versus  the  change  in  pH 
as  shown  in  Figure  9.   Superimposed  on  both  plots  is  a  hypothel  ica"1  solul 
buffered  only  by  the  aqueous  carbonate  system  with  an  alkalinity  of  106 
ueq/L  and  an  initial  pH  of  6.8,  representing  an  average  lake  for  : h '    !  d 
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The  carbonate  system,  aluminum  hydrolysis  species,  and  weak  organic 
acids  have  been  reported  as  significant  buffer  systems  for  dilute  lakes 
(Driscoll,  1980).   The  carbonate  system  is  dominant  above  pH  values  near 
5.3.   Aluminum  hydrolysis  species  and  weak  organic  acids  often  become 
dominant  buffers  below  this  point  (Driscoll,  1980). 

The  theoretical  buffer  intensities  of  each  of  these  acid-base  systems 
are  independently  presented  in  Figure  10.   To  develop  this  figure  the 
following  assumptions  were  made: 

1)  The  carbonate  system  is  in  open  equilibrium  with  the  atmosphere  al 

Pm   7  10   '   atm, 
C02 

2)  The  water  is  in  equilibrium  with  an  aluminum  trihydroxide  solid 

8.86 


phase  with  a  solubility  constant  K   =  10 


empirically  determined  for 

the  study  lakes  (see  Aluminum,  Section  4.9.1) 

3)  DOT  (average  concentration  'eve!  3.3  mg/I.)  '   nodi  1  d  as  -: 

acid  with  an  acidity  constant,  Ka  ,   =  10     •   A  similar  nodel    organic  a<  id 
has  been  empirically  verified  by  Driscoll  et  al.  (1984) 

4)  Ionic  strength  activity  corrections,  and  fluoride  and  sulfate 
complexes  with  aluminum,  are  negligible. 

From  Figure  in,  in  the  range  of  pH  values  encountered  ::;  Lhi      ;:  . 
the  dominant  buffering  systems  are  carbonate  and  water.  Aluminum  hydrolys;  ; 
species  could  provide  significant  buffering  intensity  at  low  pU  values. 
Typical  DOC  levels  will  probably  not  provide  significant  buffering. 
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5.0   DISCUSSION 

5.1   Sensitivity  Criteria 

Several  investigators  have  proposed  criteria  for  determining  lake 
sensitivity  to  acidification  based  on  pH,  alkalinity,  and  conductivity 
parameters.   Hendry  et  al .  (1980)  proposed  that  lakes  having  alkalinjties 
less  than  200  ueq/L  are  highly  sensitive,  greater  than  500  ueq 'L  alkalinity 
are  nonsensitive ,  and  between  200  to  500  ueq/L  alkalinity  are  moderately 
sensitive.  The  Ontario  Ministry  of  the  Environment  (1901)  proposed  thai 
lakes  having  less  than  40  ueq/L  alkalinity  are  extremely  sensitive  and  '■■ 
than  200  ueq/L  are  moderately  sensitive. 

Using  these  criteria  for  surface  samples,  four  lakes  (23°o)  may  be 
classified  as  extremely  sensitive,  ten  lakes  (63%)    highly  sensitive,  one 
(6%)  moderately  sensitive  and  one  lake  (Allen  T,i!^'.  G°<0  ionse   '•'■  •   PeoM 
samples  distributed  similarly  to  surface  samples  with  two  lakes  (13%)  ex- 
tremely sensitive,  nine  lakes  (G0°o)  highly  sensitive,  two  (if!0,,)  lakes 
moderately  sensitive  and  two  lakes  (Golden  and  Allen,  139o)  nonsensitive. 
Zimmerman  and  Harvey  (1979)  proposed  that  sensitive  lakes  are  characterized 
by  pH  i   6.3,  conductivity  <  40  umbos  'cm,  and  alkalinity  £  300  ueq/L.   Seven 
lakes  had  surface  and  two  lakes  had  both  surface  and  depth  samples  within 
these  criteria.  Only  one  lake  (Allen)  and  two  of  three  Golden  Lake  bottom 
samples  exceeded  all  three  criteria. 

By  all  proposed  classifications  the  study  lakes  are  preponderantly 
highly  sensitive  to  acidic  inputs.   Only  one  lake  (Allen)  is  classified  as 
not  sensitive.   Table  6  below  summarizes  lake  surface  sample  data  compared 
to  these  sensitivity  criteria. 
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Table  6 

Comparison  of  Surface  Sample  Data  to  Proposed 
Sens  i tivj  ty  Cr j  teria 


* 

Pro  po s  ed _  C r  i  t  e r i a 

Lake        Ref.  1  Ref.  2  Ref .  3 


Golden 

high 

Mowich 

high 

Unnamed 

high 

Shriner 

h  i  gh 

Crystal 

high 

Allen 

moderate 

Loui se 

h  igh 

Shadow 

high 

George 

high 

Marsh 

high 

Tipsoo 

moderate 

Frozen 

high 

Emi  ice 

high 

Reflection 

high 

Eleanor 

high 

Fan 

high 

pH   Cond  [A Ik] 


moderate 

moderate 

moderate 

modera  te 

moderate 

nonsensi  t  i  ve 

ext reme 

modera t e 

moderate 

moderate 

moderate 

extreme 

extreme 

moderate 

moderate 

extreme 


Refer  to  text  for  explanation  of  criteria. 

high  -  highly  sensitive         References: 

moderate  =  moderately  sensitive   1.  Hendry  et  al . ,  1980 

extreme  =  extremely  sensitive   2.  Ontario  Ministry  of  the  Rnvironmenl ,  19- 

-  -   exceed  (nonsensi tive)      3.  Zimmerman  and  Harvey,  1979 

-  =  below  ( sons  i  t  i ve ) 
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5.2   Temporal  and  Spatial  Trends 


Both  biogenical ly  produced  alkalinity  and  replacement  of  dilute 
snowmelt  with  ground  water  have  been  postulated  as  the  reasons  for-  increase* 
alkalinity  over  the  summer  months.   Of  the  six  lakes  temporally  studied, 
five  showed  an  increase  in  DOC,  indicating  an  increase  in  biological 
activity.   However  the  lack  of  an  alkalinity  to  DOC  correlation  indicates 
that  the  variation  in  alkalinity  is  not  highly  dependent  on  biological 
activity.   A  substantial  amount  of  the  alkalinity  increase;  in  surface 
samples  may  be  due  to  entrainment  of  higher  alkalinity  water  from  depth. 

Sulfate  ion  concentrations  increased  throughout  the  summer.   It  is  not 
possible  to  attribute  the  sulfate  increases  solely  to  either  natural 
internal  sources  or  atmospheric  deposition.   The  internal  sinks  and  sources 
of  sulfate  that  operate  on  a  seasonal  basis  tend  to  increase  sulfate 
concentrations  over-  the  summer  (Lazerte  and  Dillon  1984;  Kelly,  1984). 
Rainfall  concentrations  of  sulfate  increase  significantly  throughout  the 
summer  in  the  Washington  Cascades  due  to  tin.-  more  complete  conversion  and 
more  effective  washout  of  acid  aerosols  in  summer  rain  (Logan  -f  a  I  .  ,  1983] 
The  August  and  September  sampling  trips  followed  rainfall  even's  and  [Sf:,~ 
variation  may  in  part  reflect  rainfall  contributions.   A  change  in 
tration  with  depth  is  not  well  defined,  suggesting  that  mixing  and 
entrainment  are  not  responsible  for  sulfate  trends. 

Welch  et  al .  (1985)  proposed  using  alkalinity/cation  and 
alkalinity/conductivity  ratios  to  indicate  whether  the  lowered  spring 
alkalinity  was  simply  a  dilution  with  low  alkalinity  snowmelt  or  was  in  part 
(\ue    to  titration  of  alkalinity  with  acid  (HoS0,).   If  alkalinity  is  not 
being  replaced  then  these  ratios  will  remain  constant.   IT  replacement 
occurring  due  to  snowmelt  ibis  should  be  reflei  ted  'i  lowered  ratio:-,  in  the 


spring.   These;  ratios  remained  constant  for  the  study  lakes  and  therefore  do 
not  indicate  any  bicarbonate  replacement  with  sulfate.   The  range  of  both 
the  alkaJinity/cation  (0.27-0.33)  and  alkal inity/conductivity  (2.2  3.1) 
ratios  found  for  the  affected  Alpine  Wilderness  Area  lakes  throughout  the 
summer  are  significantly  lower  than  those  for  temporally  sampled  lakes  in 
this  study  (Alkalinity/  cation  0.92-0.91,  alkal in i ty/cond  6.6-7.9). 
Approximately  80°o  of  the  stored  acidic  ions  are  released  in  the  .firsl  20"o  of 
snowmelt  (Cadle  et  al . ,  1985).   The  resulting  spike  of  acidification  has 
been  noted  by  many  researchers  (Klock  and  Lefohn,  1935;  Welch  et  al . ,  !935). 
Since  lake  sampling  in  this  study  began  after-  the  initial  spring  snowmelt, 
the  potential  acid  spike  was  probably  not  included  in  the  sample  data  and  it 
is  therefore  questionable  whether  the  effect  would  have  been  detected. 

5.3  Buffer  System 

In  lakes  dominated  by  the  carbonate  system,  resistance  to  pH  depression 

is  achieved  through  the  loss  of  alkalinity.  For  the  range  of  pH  values 

encountered  in  this  study  the  maximum  buffering  capacity,  as  shown  by  the 

peak  in  Figure  9,  occurs  at  pH  =  pK   of  the  carbonate  system  near  pH  = 

a  I 

6.4.   The  minimum  buffering  capacity  occurs  at.  the  ro^  equivalence  point 
near  pH  =  5.25  for  dilute  systems.   Therefore,  systems  at  a  pH  near  pKa 
would  require  relatively  much  higher  acid  input,  or  would  exhibit  a  much 
greater  change  in  alkalinity,  for  a  given  change  in  pH  than  systems  nearer 
the  CO   equivalence  point. 

The  similarity  between  the  lake  data  and  the  theoretical  curve  with 
regard  to  equivalence  points  and  slopes  indicates  the  dominance  of  the 
carbonate  buffer  system  in  the  lakes  sampled.   The  pKa  value,  or  poin   >'" 
maximum  buffer  intensity,  is  sensitive  to  changes  in  ion.ii     ength  ■•  ' 
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temperature.   Although  the  pH  value  of  the  pKa   point  docs  nol  change  as  ^ 

function  of  carbonate  concentration,  the  magnitude  of  the  buffering 

intensity  does.   In  contrast,  the  buffer  intensity  minimum  defined  by  the 

C02  equivalence  point  varies  both  in  magnitude  and  pH  value  with  changes  in 

total  carbonate  concentration,  or  p       This  fact  may  explain  the  better 

theoretical  fit  at  the  maximum  than  at  the  minimum  buffer  intensity  points. 

Lake  pH  values  tend  to  stabilize  where  a  buffering  capacity  is  high  rather 

than  at  minimum  points,  creating  a  bimodal  distribution  (Wright  and 

Gjessing,  1976;  Zimmerman,  1982).   Zimmerman  (1982)  proposed  that  this 

distribution  may  be  used  to  indicate  acidification.   Nonacidified  areas 

would  have  lake  pH  values  distributed  around  pK  ,  (pH  =  6.4).   Acidified 

a  l 

areas  would  show  an  increase  in  lake  pH  distribution  where  natural  water  or 
aluminum  buffering  potential  increases,  near  pH~5 .   Lakes  in  the  transition 
zone  ( pH  =  5.25-5.7)  are  characterized  by  severe  pH  fluctuations  and  are 
extremely  sensitive  to  seasonal  pulses  of  acid,  as  during  snowmeU 


distribution  due  to  acidification  is  noticeable.   Lakes  would  remain  in  the 
transition  zone  until  almost  all  of  the  alkalinity  is  consumed  (Henricksen, 
1980).   The  acidified  areas  of  the  Adirondacks,  when  compared  !.o  h.Lsl'u  \  ' 
data,  have  shown  a  definite  shift  to  the  acidified  region  of  the  bimodal 
distribution  (Driscoll,  1985).   The  lakes  in  this  study  showed  a  pH 
distribution  around  the  maximum  buffer  intensity  (pK   )  point  (median  pil  = 
6.4,  Table  4),  or  nonacidified  region. 

The  apparent  lack  of  an  additional  weak  buffering  system  indicates  that 
equilibration  at  a  higher  p„n   is  the  best  explanation  for  values  falling 

o 

above  the  theoretical  pH  vs.  alkalinity  line  in  Figures  3  and  '.    The  average 

•3.01 


no 


as  calculated  from  lab-equilibrated  pH  and  alkalinity,  was  in 


50 


atm.   This  value  is  3.1  times  higher  than  the  global  average  p    at  sea 

-3  5  - 

level  of  10  '  '   atm.   Similar  results  were  reported  for  the  Central 

Washington  Cascade  lakes  and  Oregon  Cascade  lakes  (Logan  el  al .  ,  1983; 
Nelson  and  Delwiche,  1933).   Logan  et  a"!.,  (1983)  attributed  the 
equilibration  at  a  higher  p    to  temperature  differences  between  lakes 
sampled  and  the  laboratory.   Nelson  and  Delwiche  (1983)  argued  that  since 
pH  stabilizes  asymptotically  over  long  periods  of  time,  the  existence  of  a 
kinetic  limitation  on  the  equilibration  of  the  water  sample  with 
atmospheric  p~n   is  indicated.   Such  a    limitation  may  explain  the  apparent 
supersaturation  caused  by  the  higher  soil  p    relative  to  thai  of  the 
atmosphere  (Elzerman,  1983).   The  supersaturation  illustrated  in  Figure  -1 
probably  reflects  this  effect. 

5.4   Aluminum  Chemistry 

The  total  aluminum  concentrations  determined  in  this  study  were  similaj 
to  the  total  aluminum  concentrations  found  by  the  U.S.G.S.  (1986)  study 
Elevated  levels  of  total  aluminum  are  influenced  by  both  organic  and 
inorganic  complexation  and  are  not  a  response  to  acidification.   Seasonal 
variation  in  total  aluminum  is  a  response  to  varying  levels  of  organic  ill} 
complexed  aluminum. 

The  general  pattern  of  aluminum  spec.iat.ion  observed  in  the  study  iak*  s 
agrees  well  with  theoretical  thermodynamic  equilibrium  models.   The  lake 
water  data  resembles  aluminum  trihydroxide  equlilibria  and  indicates  thai 
the  relatively  higher  levels  of  aluminum  found  in  some  lakes  is  not 
associated  with  higher  free  ionic  aluminum  ([Al°  ])  concenlrul  ions  or 
acidification.  Vor   the  lakes  studied,  nearly  all  of  the  labile  ill  :ti  \  ■ 
theoretically  predicted  to  be  in  the  A1(0H)    form  sine1  the  pll   a  '  ■       wore 
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generally  above  6.  Temporal  changes  in  total  aluminum  concentration  appear 
to  be  a  function  of  changes  in  the  organically  complexed  aluminum  component 
which  responded  to  changes  in  DOC. 

Although  dissolved  organic  carbon  levels  are  low  in  the  study  lakes 
(mean  DOC  =  1.3  mg/L),  the  average  organic  aluminum/total  aluminum  ratio  of 
0.29  indicates  that  a  significant  amount  of  the  aluminum  is  organically 
bound.   The  nonlabile  or  organically  complexed  aluminum  is  strongly 
correlated  to  DOC  as  seen  by  equation  4  2  (App.  2).   Simitar  relationships 
have  been  found  by  other  researchers  (Wright,  1985;  Driscoll,  1985;  a:)(] 
Johnson,  1981).   Levels  of  organic  aluminum  complexation  are  influenced  bj 
the  altitude  and  position  of  the  lakes  in  the  watershed,  the  upper  reaches 
containing  less  organically  complexed  aluminum. 

Both  the  early  summer  and  fall  sampling  trips  occurred  during  runoff 
events  of  snowmelt  or  rainfall.   Such  events  increase  the  levels  of  DOC 
organically  complexed  aluminum  in  streamflows  (Driscoll,  1984;  Johnson  et. 
al . ,  1934).   It  is  likely  that  increased  levels  of  aluminum  in  September  are 
a  response  to  inflow  rather  than  to  biogenically  derived  DOC.   The  one  ex- 
treme total  aluminum  value  of  80  ug/L  for  Tipsoo  Lake  is  much  lower  than  the 
average  of  382  ug/L  for  three  acidified  Adirondack  lakes  (Driscoll,  !983). 
The  mean  total  aluminum  lake  surface  concentration  of  11  ug/L  is  similar  to 
the  mean  of  10  ug/L  cited  as  background  levels  for  natural  water  and  less 
than  the  16  ug/L  for  Oregon  Cascade  lakes  (Stumm  and  Morgan,  1981;  Nelson 
and  Delwiche,  1983) . 

3.5   Sulfate 

To  relate  sulfate  concentrations  to  redui  I  ion  in  alkalinity  or  pU 
depression,  Wright  (1983)  developed  an  "acidity  index"  equation  based  on 
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theoretical  ion  balance  equations: 


([H'J    [HC0o  ]) 


-0.91  ([Ca2:]  -  [Mg2 


[SO,2] 


An  increase  in  [SO    ]  through  acid  deposition  causes  either  an  increase  in 

4 

( [H  ]  -  [HCO   ]),  which  is  equivalent  to  a  decrease  in  alkalinity,  an 
increase  in  ([Ca   ]  +  [Mg'"  ]),  or  both. 

Not  all  sulfate  can  be  attributed  to  anthropogenic  deposition. 
Concentrations  of  background  sulfate,  attributed  in  part  to  terrestrial 
dust,  exist  in  pristine  lakes  (Wright,  1983).   Background  sulfate  !  eve'l  s  "r 
20  ueq  'L  and  15  ueq/L  proposed  by  Galloway  et  al.  (1983)  and  Wright  •.  ■' 
Henricksen  (1983),  respectively,  and  estimates  from  Wright's  background 
(Bkg)  sulfate  equation; 

3kg  [SO,2"]  =  0.09  ([Ca2  ]-[Mg2"])-[X^]-[K^ 

are  greater  than  those  typically  encountered  in  this  study  X  =  00 . 4  ueq. 
Estimates  from  an  equation  developed  for  this  study  using  Wright's  procedure 


Bkg  [SO,"  ]  -  7.1-0.05  E  [cations] 

4 


were  similar  to  observed  total  [SO    ]  (mean    10.00  ueq/1,  median    12 

4 

ueq.'L)  suggesting  levels  at  or  near  background  levels. 

Wright  (1983)  found  net  [SO  2~]  (total  [SO  2~]    Bkg  [SO  "])  to  be 

closely  correlated  (r2  =  0.99)  with  alkalinity  loss  (  A[AJk])  for  North 

American  lakes.   Hendry  and  Brezonik  (1983)  and  Zimmerman  (1902),  using 

similar  indices  (Net  [SO,2']  +  [Alk],  Net  [SO.  "]/  [Alk]),  found  a  strong 

4  4 

negative  correlation  with  pH  (r   =  0.71).   These  methods  did  not  associato 
net  [SO  ~    j  with  loss  of  alkalinity  or  oil  depression  in  the  lakes  shi  !.ied 
A  number  '^f    researchers  have  documented  the  mitigating  efforts  of 
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natural  processes  such  as  soil  leaching  and  biochemical  reduction  on  sulfate 
deposition.   Correlating  sulfate  deposition  rates  to  aquatic  concentrations 
without  an  extensive  mass  balance  analysis  is,  therefore,  somewhal  tenuous. 
Sulfate  ion  flux  is  increased  by  acidic  deposition,  resulting  in  an 
alkalinity  decrease,  and  is  consistent  with  the  titration  concept  (Lazerte 
and  Dillon,  1984) . 

The  median  lake  surface  [SO  ~    for  the  lakes  studied  of  8.3  ueq/L  is 

4 

at  the  low  end  of  the  8.3  ueq/L  to  14.8  ueq/L  range  in  precipitation 
reported  for  Snoqualmie  Pass,  Washington  (Logan  et  al  .  ,  1983).   When  high 
local  sources  of  sulfate  (Allen  Lake)  are  nut  considered  in  this  study,  the 
average  concentration  (8.4  ueq/L)  is  greater  than  the  average  for  the  Oregon 
Cascade  lakes  (5.8  ueq/L)  and  less  than  the  average  for  the  Alpine  Lakes 
Wilderness  Area,  Washington  (13.]  ueq/L)  (Nelson  and  De!wich<  .  1988 
1933).   Assuming  similar  base  geology,  the  disparity  in  sulfate 
concentrations  indicates  that  some  sulfate  in  Washington  is 
anthropogenically  derived  and  that  the  effect  is  less  important  at  Mount 
Rainier  than  in  the  Alpine  Wilderness  Area  farther  north.  The  level  of 
acidification  indicated  by  the  disparity  in  sulfate  data  is  below  the 
resolution  of  the  models  used  in  this  study  to  determine  the  extei  I  uf 
acidification  (see  section  5.6).   The  effect  of  sulfate  deposition  on  lake 
water  chemistry  would  be  greatest  during  the  initial  snowmelt  and  could 
result  in  significant  peaks  of  acidification  in  the  extremely  sensitive 
lakes  (Table  5).   Further  study  would  be  required  to  determine  the 
presence  and  extent  of  the  possible  snowmelt  induced  pulse  of 
acidification.   The  trend  of  increasing  anthropogenic  sulfate  going 
north  is  consistent  with  prevailing  weather  patterns  foi  Washington 
State. 
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5.6   Model  Applications 
5.6.1  Electroneutrality 

To  determine  the  extent  of  acidification  in  the  study  lakes,  water 
chemistry  data  is  analyzed  using  models  which  are  developed  from  equations 
presented  in  Table  1.   The  theoretical  and  empirical  equations  from  Table  1 
also  provide  a  basis  for  comparing  this  studies  data  to  that  of  other 
researchers . 

By  assuming  carbonic  acid  weathering  and  fulfilling  the  requirements  <.r 
electroneutrality,  the  correlation  between  alkalinity  and  total  ca! 
should  result  in  a  slope  =  1  and  y  intercept  •  0  ("<;.  7,  Table  1). 
strong  mineral  acid  weathering  has  occurred  in  a  given  area,  the  slope  will 
be  less  than  one,  indicating  the  release  of  cations  without  accompanying 
alkalinity,  or  the  titration  of  bicarbonate   The  closeness  to  theo-y  >* 
alkalinity  vs.  cation  equation  (#34;  [Alk]  =  -7.7  -  0.05  :<    [cal  ions]   \v 
2)  indicates  that  carbonic  acid  weathering  dominates  the  lake  systems 
studied  and  that  little,  if  any,  acidification  has  occurred.   A  very  similar 
equation  ([Alk]  =  -8.66  +  1.05  x  [cations])  was  determined  for  the 
nonacidified  Oregon  Cascade  lakes  (Nelson  and  Delwiche,  1983). 

Burns  (1981)  proposed  an  empirical  method,  based  on  the  theoretical Ij 
derived  electroneutrality  equation,  by  showing  that  in  areas  unaffected  by 
acid  deposition,  the  slope  of  the  carbonic  acid  weathering  line  ([Alk]  vs. 
[cations])  was  greater  than  or  equal  to  0.75.   A  slope  of  one  is  seldom 
attained,  which  is  attributed  to  either  an  alternative  source  of  weak  acid 
or  the  omission  of  a  component  anion  (Burns,  1981;  Zimmerman,  1982).   The 
ratio  for  acidified  areas  was  significantly  lower,  i.e..  0.064,  Scandanavi;  . 
0  185,  Hubbard  Brook,  NT.K.;  and  0 .  562 .  North  Carolina.   Figin>e  11 
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illustrates  Burns'  model  with  the  study  lakes  superimposed.   The  study 
lakes  could  lose  an  average  of  11";  of  their  bicarbonate  alkalinity  which  is 
equal  to  an  average  of  13.7  ueq/L  of  acidification  for  surface  samples,  and 
not  appear  acidified  according  to  Burns'  method. 

Only  one  sample,  (Frozen  lake,  depth,  August)  falls  below  the  empirical 
[Alk]  vs.  ^[cations]  ratio  0.75  defined  by  Burns  as  the  boundary  between 
natural  conditions  and  acidification.   This  error  may  be  due  to  the  extreme 
diluteness  of  Frozen  Lake  (Zimmerman,  1982).   However  the  ratio  for  Frozen 
lake  surface  sample  lies  above  the  limit.   Since  the  lake  was  not 
stratified,  the  averaged  values  for  the  parameters  were  used  which  1  :  <:  ibovo 
the  cutoff  (0.75).   The  low  dissolved  silica  and  aluminum  concentrations 
the  presence  of  a  snowpack  near  the  lake  when  sampled  in  Angus!  indicate 
that  Frozen  Lake  receives  most  of  its  water  from  precipitation  runoff  rather 
than  groundwater.   The  low  pH  (5.7)  and  extremely  sensitive   iture 
(conductivity  =  3.5  umho/cm,  alkalinity  =  15.5  ueq/L)  are  within  the 
expected  range  for  precipitation  dominated  lakes  (Klok  and  i.efohn,  1935). 
Therefore,  by  Burns  method,  none  of  the  sampled  lakes  has  received 
significant  acid  inputs. 

To  empirically  determine  if  the  key  lakes  "elevated"  aluminum  could  be 
associated  with  low  levels  of  acidification  [Alk]  vs .  E  [cation]  equations 
were  developed  for  both  high  aluminum  (Golden,  Mowich,  Shriver)  and  typical 
aluminum  (Mowich,  Crystal,  Allen)  lakes.   A  higher  level  of  acidification 
would  result  in  a  lowered  intercept  and  slope  in  the  [alk]  vs.  Z[cation] 
correlation.   The  high  aluminum  lakes  correlation's  negative  intercept  and 
large  slope  ( eq .  39)  appear  to  be  an  artifact  of  the  much  higher  ionic 
strength  of  the  Golden  Lake  depth  samples.   When  Golden  Lake  depth  samples 
ire  omitted,  the  high  aluminum  (Eq.  }Ci)    and  typical  aluminum  (Fo.  H) 
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equations  are  not  significantly  different  (p  <:  0.03).   Thus  relatively  high 
aluminum  levels  in  the  high  aluminum  lakes  do  not  appear  to  be  the  result 
of  acidification. 

As  discussed  earlier  the  y-intercept  of  the  electroneutral Aty  equations 
(Eqs.  6  anni trogen-oxide-derived  acidification 
(nitric  acid)  is  unlikely  to  have  occurred  in  the  study  area. 

Models  based  on  [Alk]  vs.  [Ca2+]  (Henricksen,  1979)  and  [Alk]  vs.  [Ca2^] 

2  +■ 
+  [Mg   ]  (Henricksen,  1980)  are  essentially  empirical  simplifications  of  the 

theoretical  electroneutrality  equation  (Fqs.  10.  11.  Table  1).   The 

empirical  relationships  derived  from  unacidified  Norwegian  lakes  '  .-,   de  i 

original  alkalinity  underestimates  observed  nitrogen-oxide-derived  acidif icatlo; 

(nitric  acid)  is  unlikely  to  have  occurred  in  the  study  area. 

Models  based  on  [Alk]  vs.  [Ca2"]  (Henricksen,  1979)  and  [Alk]  vs.  [Ca2~] 

~  [Mg""  ]  (Henricksen,  1980)  are  essentially  empirical  simp] i flea '  '  >.:.  of  '  '.;< 

theoretical  electroneutrality  equation  (Fqs.  10,  11.  Table  1).   The 

empirical  relationships  derived  from  unacidified  Norwegian  lakes  to  describe 

original  alkalinity  underestimates  observed  alkalinity  in  the  Mt .  Rainier 

2  -        2- 
study  lakes  by  an  average  of  17%  using  the  alkalinity  vs.  ( [Ca   ]  -  (Mg-  ] 

o ... 

model  and  23 . 8%   using  the  alkalinity  vs.  [Ca-  ]  model.   Multiple  regression 
indicates  that  this  is  due  to  the  significant  alkalinity  associated  with 
potassium  and  sodium.   Welch  and  Chamberlain  (1982)  found  that.  Henricksen 's 
(1979)  alkalinity  vs.  [Ca   ]  model  underestimated  observed  alkalinity  by  22°« 
for  Washington  lakes  in  general. 

3.6.2   Ion  Concentration  vs.  pH  Plots 

2 
Henricksen  (1979)  proposed  a  pH  vs.  [Ca   ]  equation  to  distinguish 


between  acidified  and  nonacidified  waters  (e.g.  20, 


e  l ) .  i he  bourn 
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line  ([Ca   ]  vs.  pH )  between  acidified  and  nonacidi  f  led  lakes  was  empiri 
cally  determined  using  data  from  59  nonacidified  Norwegian  lakes 
(Henricksen,  1979).   Figures  12  and  13  present  Pa  vs.  pH  data  far  the  study 
lakes  using  lab-equilibrated  and  field  pH  values,  respectively,  and  with 
Henricksen's  boundary  superimposed.   Some  loss  of  alkalinity,  or  a  moderate 
drop  in  pH ,  could  occur  from  acidification  and  most  lakes  would  appear 
unimpacted  according  to  Henricksen's  model.   Figure  14  compares  the  study 
lakes  with  the  Washington  Cascade  lakes  studied  by  Logan  (1982). 

Comparison  of  Figures  12  and  13  shows  that  selecting  field  pK  over  1  sb 
equilibrated  pH  greatly  increases  the  apparent  level  of  acidification  f. 
the  lakes  studied.   Being  an  intensity  parameter,  ;>T! ,  when  nieasu   '  in  ".■ 
field,  fluctuates  significantly  due  to  p    variation,  photosynthesis, 
respiration,  wave  action,  temperature,  and  measuring  method.   These  Factors 
are  mr'nimized  by  determining  pH  under  consistent  !abo  itorj  ■   •  " '  ' 
However,  since  no  defined  standard  exists,  it  is  difficult  to  compat   p." 
data  among  investigators. 

Kramer  and  Tessier  (1982)  suggest  that  plots  of  total  cation 
concentration  vs.  pH  provide  a  more  suitable  model  for  testing  acidification 
since  pll  and  alkalinity  are  directly  related  and  rh«  boundary  l\\       has  u!.< 
retical  justification,  (Eq.  10,  Table  1).   Similar  to  previous  methods  this 
method  indicated  no  acidification  in  the  study  lakes. 
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5.J1     CONCLUSIONS 

With  respect  to  the  stated  objectives,  the  following  conclusions  may  be 
drawn  from  this  study: 

1)  Mount  Rainier  National  Park  lakes  are  of  a  generally  dilute  ionic 
composition  relative  to  average  North  American  freshwater  lakes  but  are 
typical  of  Cascade  Mountain  Lakes. 

2)  Most  lakes  in  Mount  Rainier  National  Park  are  highly  susceptible 
to  acidification  from  acid  precipitation  due  to  their  diluteness  and  poor 
buffering  capacity. 

3)  Observed  dissolved  aluminum  concentrations  simulate  theoretical 
equilibrium  with  an  aluminum  trihydroxide  phase  similar  to  natural  gibbsite. 
The  relatively  higher  concentrations  of  total  aluminum  observed    so>n 
lakes  is  not  associated  with  the  free  aluminum  ion  cited  as  a  symp:  ■• 
acidification,  but  rather  is  attributable  to  natural  soluble  complexes  with 
inorganic  and  organic  ligands. 

4)  Based  on  ionic  composition,  Mount  Rainier  lakes  exhibit  no 
detectable  level  of  acidification. 
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Empirical  Correlations  of  Selected  Chemical  Parameters  Using  Averaged  Parameter 

Values,  Sixteen  Study  Lakes,  Mount  Rainier  National  Park  (Units  are  mole/1  or 
equivalent/1  unless  otherwise  noted.) 
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Equation  (b  +  mx)  r_ 

0.0  +  1.08  pH  Field  0.54 

-3.07  +  0.74  pH  0.81 

-3.9  +  0.69  pH  0.74 

-1.94  +  0.22  pH  0.31 

-2.73  +  0.46  pH  0.22 

-2.12  +  0.62  pH  0.77 

-3.01  +  0.73  pH  0.80 

-2.46  +  0.50  pH  0.26 

-1.56  +  1.10  pH  0.06 

0.24  -  0.14  pH  0.25 

-2.52  +  0.67  pH  0.78 

0.69  +  0. 18  pH  Field  0.21 

-2.73  +  0.33  pH  0.22 

-3.27  +  0.45  pH  0.20 

-2.54  +  0.55  pH  0.71 

0.005  exp  (1 .19  pH)  0.73 

-100  +  17.9  pH  0. 17 

-148  +  25.1  pH  Field  0.31 

-2.75  +  0.8  pH  0.01 

-25.5  +  7.43  Cond  0  99 

-2.04  -  0.52  Cond  0.93 

9.77  +  0.61  Cond  0.67 

-0.22  -  0.11  Cond  0.43 
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-7.71 


0.95  [Cations] 


184    +   22.75  [SiO  1 

2 J mg/1 


104.5  +   34.80  [Al] 


ug/1 


105.4  +   0.54  DOC 


mg/1 


-9.1    +   1.07  [Cations]-2[S0.   ] 

4 

-20.88  +   1 .05  [Cations] 

-0.5    +   0.85  [Cations] 
depth  sample) 

13.01  +   0.82  [Cations] 


2.58 


4. 15  DOC 


mg/1 


2  + 


-29    +   1.32  [Ca^  ] 
(not  seasalt  corrected) 
10    +   1.23  [On   ] 


0.0    +   1.21  [CI 


c  + 


0.99 
0.46 
0.21 
0.58 
0.99 
0.29 
0.57 
0.98 

0.97 
0.92 
0.73 
0.53 
0.97 
0.99 

0  .  59 
0.35 
0.53 
0.99 

1  .00 
0.95 
1  .00 
0.80 


•14     +   0.91  ([Ca2+]+[Mg2+]) 


73 


Appendix  3:   Data  Summary  for  Mt .  Rainier  Study 

Table  A3-1:   Surface  Samples 

Aluminum,  ug/1 
PH       . 


Lake    Trip  Field   Lab   Total    Organic   Labile 


Golden 

1 

5.42 

6.27 

35 

96 

5.97 

29.99 

2 

6.48 

6.61 

21 

17 

4.23 

16.93 

3 

6.55 

6.77 

42 

99 

17.90 

25.09 

Mowich 

1 

6.86 

7.05 

2 

74 

0.41 

2.33 

2 

6.20 

7.09 

1 

42 

0.33 

1  .09 

3 

6.56 

7.08 

6 

39 

0.26 

6.13 

Unnamed 

1 

5.85 

6.30 

1 

14 

0.51 

0.63 

2 

5.78 

6.38 

4 

84 

0.61 

4.23 

3 

6.27 

6.75 

4 

83 

2.  10 

2.73 

Shriner 

1 

6.06 

7.09 

14 

00 

0.41 

13.59 

2 

6.28 

7.18 

16 

51 

0.32 

16.19 

3 

6  .  35 

6.78 

16 

58 

1.47 

15.11 

Crystal 

1 

5.96 

6.59 

1 

57 

0.30 

1.27 

2 

6.60 

7.13 

4 

60 

0.24 

4.35 

3 

6.55 

6.65 

9 

38 

5.57 

3.80 

Allen 

1 

6.60 

7.58 

2 

37 

1  .23 

1.08 

2 

6.80 

7.59 

3 

05 

0.54 

2.51 

O 

o 

6.71 

7  .  55 

4 

94 

3.46 

1  .48 

Louise 

1 

6.50 

6.37 

1 

73 

0.50 

1  .23 

Shadow 

1 

6.50 

6.71 

2 

81 

1  .48 

1  .33 

George 

1 

7.01 

7.  15 

2 

26 

0.37 

1.89 

Marsh 

2 

6.31 

6.14 

2 

83 

1.38 

1.45 

Tipsoo 

2 

7.15 

7.03 

75 

03 

5.93 

69.  10 

3 

7.43 

6.76 

81 

03 

10.77 

70.26 

Frozen 

2 

5.33 

5.73 

0 

97 

0.38 

0.59 

Eunice 

2 

5.71 

6.35 

4 

00 

1  .  10 

2.90 

Ref lectio 

3 

6.20 

6.73 

6 

41 

4.01 

2.40 

Eleanor 

O 

6.50 

6.78 

5 

15 

1  .24 

3.91 

Fan 

X  'A 

5.62 

7 

67 

4.01 

3.66 

max 

7.43 

7.59 

81 

03 

17.90 

70.26 

min 

0.00 

5.62 

0 

97 

0.24 

0.59 

avg 

6.  16 

6.75 

13 

25 

2.66 

10.60 

std 

1.25 

0.48 

20 

21 

3.77 

17.73 

ANC 


Anions,    mg/1 


ueq/1 

mg/] 
as 

2- 

-- 

- 

2  - 

- 

Lake      T 

rip 

CaCO 

co 

HC0 

CI 

SO 

F(ug/ 

3 

3 

3 

4 

Golden 

1 

57 

57 

2.88 

0 

00016 

3.54 

0.46 

0.37 

11 

20 

2 

58 

06 

2.90 

0 

00036 

3.56 

0.46 

0.50 

11 

11 

3 

61 

20 

3.06 

0 

00054 

3.74 

0.54 

0.61 

9 

47 

Mowich 

1 

79 

16 

3.96 

0 

00134 

4.83 

0.41 

0  .  35 

11 

20 

2 

84 

53 

4.23 

0 

00157 

5.  16 

0.55 

0.39 

90 

3 

31 

81 

4.09 

0 

00148 

5.00 

0.49 

0.40 

9 

40 

Unnamed 

1 

55 

52 

2.78 

0 

00017 

3.42 

0.46 

0.50 

13 

27 

2 

66 

24 

3.31 

0 

00024 

4.07 

0.43 

0  .  57 

14 

96 

3 

80 

97 

4.05 

0 

00069 

4.95 

0  .  5 1 

0.65 

10 

85 

Shriner 

1 

120 

67 

6.03 

0 

00223 

7.37 

0.27 

0  .  09 

'  5 

75 

2 

132 

10 

6.61 

0 

00301 

8.06 

0.33 

0.97 

m 

30 

3 

122 

15 

6.11 

0 

00111 

7.46 

0.42 

0  .  95 

•6 

16 

Crystal 

1 

117 

35 

5.87 

0 

00069 

7.  17 

0.  10 

0.  18 

13 

27 

2 

126 

25 

6.31 

0 

00256 

7.71 

0.07 

0.20 

14 

33 

121 

30 

6.07 

0 

00082 

7.41  • 

0.  17 

0.25 

12 

^  -f 

Allen 

1 

336 

73 

19.34 

0 

02211 

23.59 

0.54 

3.23 

18 

72 

2 

418 

06 

20 .  no 

0 

02445 

25.50 

0  .  54 

5.03 

18 

08 

3 

411 

71 

20.59 

0 

02100 

25.  11 

0  .  50 

0.01 

1 8 

7 •> 

Louise 

1 

28 

91 

1  .45 

0 

00010 

1  .79 

0.35 

0.  10 

3 

71 

Shadow 

1 

49 

32 

2.47 

0 

00038 

3.02 

0.  10 

0.15 

17 

75 

George 

1 

198 

07 

9.90 

0 

00421 

12.09 

0.62 

0.31 

13 

35 

Marsh 

2 

53 

06 

2.65 

0 

00011 

3.23 

0.31 

0.79 

19 

72 

Tipsoo 

2 

234 

72 

11.74 

0 

00373 

14.32 

0.  15 

0  .  58 

23 

90 

3 

228 

44 

11.42 

0 

00198 

13.94 

0.17 

0.50 

17 

57 

Frozen 

2 

17 

45 

0.87 

0 

00002 

1.18 

0.01 

0.43 

10 

04 

Eunice 

2 

38 

14 

1.91 

0 

00013 

2.35 

0.33 

0.20 

6 

52 

Reflection 

57 

91 

2.90 

0 

00047 

3.54 

0.27 

0.42 

11 

20 

Eleanor 

\J 

1  oo 

1  oo 

72 

6.69 

0 

00121 

3  .  1 7 

0.50 

1  .37 

20 

27 

Fan 

3 

31 

25 

1  .  56 

0 

00002 

2.05 

0.32 

0.7  1 

11 

62 

max 

418 

06 

20.90 

0 

02445 

25  .  50 

0.62 

6.01 

23 

90 

min 

17 

45 

0.87 

0 

00002 

1.18 

0.01 

0.  10 

6 

52 

avg 

125 

94 

6.30 

0 

00338 

7.70 

0.36 

0.99 

13 

97 

std 

109 

28 

5.46 

0 

00671 

6.65 

0.  17 

1  .45 

4 

12 

Cations,  mg 

'] 

2  + 

~~2  + 

+ 

- 

Lake     Tr 

ip 

Ca 

Mg 

Na 

K 

SiO 

2 

Golden 

1 

0.85 

0.  12 

0.61 

0.06 

4.  16 

2 

0.91 

0.15 

0.64 

0.06 

4.10 

3 

0.98 

0.  18 

0.74 

0.07 

4.25 

iMowich 

1 

1.43 

0.16 

0.48 

0.03 

3.00 

2 

1.52 

0.  17 

0.47 

0.04 

2.93 

3 

1.48 

0.  16 

0.49 

0.04 

3.07 

Unnamed 

1 

0.80 

0.  14 

0.81 

0.07 

5.61 

2 

1  .02 

0.  17 

0.79 

0.05 

5  .  59 

3 

1  .  15 

0.  18 

0.84 

0  .  04 

5.92 

Shriner 

1 

2.16 

0.26 

0.79 

0.22 

3.51 

2 

2.50 

0.31 

0.62 

0.09 

3.26 

3 

2.22 

0.29 

0.62 

0.05 

2.25 

Crystal 

1 

1.93 

0.  14 

0.48 

0.  16 

3.90 

2 

2.25 

0.  14 

0.42 

0.06 

3.47 

3 

2.26 

0.14 

0.38 

0.07 

3.06 

Allen 

1 

7.73 

0.76 

0.94 

0.05 

5.  70 

2 

8.59 

0.86 

0.93 

0.09 

4.23 

3 

7  .  99 

0.83 

0  .  94 

0.04 

3.03 

Louise 

1 

0.47 

0.06 

0.42 

0.03 

2.43 

Shadow 

1 

0.65 

0.  11 

0.32 

0.06 

3.64 

George 

1 

3.36 

0.32 

0.85 

0.  13 

4.00 

Marsh 

2 

0.67 

0.  15 

0.58 

0.  11 

4.07 

Tipsoo 

2 

4.67 

0.23 

0.61 

0.09 

2.35 

3 

4.17 

0.22 

0.67 

0.09 

0.80 

Frozen 

2 

0.20 

0.02 

0.24 

0.02 

3.29 

Eunice 

2 

0.69 

0.07 

0.41 

0.03 

1  .92 

Reflection 

0.84 

0.  16 

0.53 

0.03 

3.30 

Eleanor 

3 

2.28 

0.43 

0.73 

0.03 

4  .  58 

Fan 

3 

0.44 

0.09 

0.47 

0.08 

3.  15 

max 

8.59 

0.86 

0.94 

0.22 

5.92 

min 

0.20 

0.02 

0.24 

0.02 

0.80 

avg 

2.28 

0.24 

0.61 

0.07 

3.62 

std 

2.25 

0.21 

0.  19 

0.04 

1  .15 

Ca  1  c  . 

Total 

Total 

Ratio 

DOC 

Cond . 

Ionic 

Cations 

Anions 

Tot. Cat 

Lake     Tr 

ip 

mg/1 

umho/cm 

Str. 

ueq/1 

ueq/1 

/An . 

mole/1 

Golden 

1 

2.30 

15 

00 

0 

00011 

79.87 

78.96 

1.01 

2 

1.  15 

10 

51 

0 

00012 

86.74 

82.01 

1.06 

3 

3.95 

11 

05 

0 

00013 

96.39 

89.46 

1  .08 

Mowich 

1 

1  .22 

14 

90 

0 

00015 

105.26 

98.56 

1  .07 

2 

0.54 

12 

10 

0 

00016 

110.49 

108.49 

1  .02 

3 

1.21 

11 

00 

0 

00015 

108.11 

104.51 

1  .03 

Unnamed 

1 

0.53 

12 

90 

0 

00012 

87.54 

79.72 

1  .  10 

2 

0.54 

12 

10 

0 

00013 

98.97 

91  .09 

1  .09 

O 

0.87 

10 

60 

0 

00015 

108.83 

109.52 

0  .  99 

Shriner 

1 

0.67 

19 

10 

0 

00023 

167.07 

14  3.37 

1.17 

2 

1.32 

16 

80 

0 

00026 

177.36 

162.44 

1  .  09 

3 

1  .31 

18 

90 

0 

00024 

160.89 

154.50 

1  .04 

Crystal 

1 

1  .36 

17 

30 

0 

00018 

130.36 

124.65 

1  .05 

2 

1  .24 

18 

05 

0 

00020 

141.53 

133. 18 

1.06 

3 

1  .28 

15 

00 

0 

00020 

140.67 

132.04 

1  .07 

Allen 

1 

0.67 

56 

00 

0 

00074 

483.53 

471.22 

1  .03 

2 

0.70 

57 

10 

0 

00086 

534.24 

551 .26 

0.97 

3 

1.25 

55 

90 

0 

00033 

501 .87 

553.69 

0.91 

Louise 

1 

0.35 

5 

95 

0 

00006 

47.03 

41  .24 

1.14 

Shadow 

1 

0.96 

7 

10 

0 

00008 

56.  17 

56.  12 

1  .00 

George 

1 

0.97 

28 

20 

0 

00033 

231.51 

222.68 

1.04 

Marsh 

2 

0.57 

10 

20 

0 

00011 

72.50 

79.15 

0.92 

Tipsoo 

2 

2.05 

31 

99 

0 

00040 

276.17 

252.23 

1  .09 

3 

3.30 

31 

00 

0 

00037 

253.39 

244.66 

1  .04 

Frozen 

2 

0.49 

O 

66 

0 

00004 

22.41 

28.32 

0.79 

Eunice 

2 

0.85 

9 

50 

0 

00008 

58.06 

51.96 

1.  12 

Reflection 

2.25 

7 

50 

0 

00011 

73.40 

74.79 

1  .  05 

Eleanor 

0.73 

13 

38 

0 

00028 

179.78 

187.75 

0.96 

Fan 

o 

3.09 

9 

40 

0 

00008 

51.21 

55.30 

0.92 

max 

3.95 

57 

10 

0 

00086 

534.24 

553.69 

1.17 

min 

0.49 

3 

66 

0 

00004 

22.41 

28.32 

0.79 

avg 

1.31 

18 

87 

0 

00024 

160.23 

157.35 

1  .03 

std 

0.88 

14 

36 

0 

00021 

132.08 

137.42 

0.08 

Ratio    Ratio   Ratio 
An /Cat  ANC/Cat  ANC/Cond 
Lake     Trip 


Golden 

1 

0.99 

0.72 

3.84 

2 

0.95 

0.67 

5.52 

3 

0.93 

0.63 

5.54 

Mowich 

1 

0.94 

0.75 

5.31 

2 

0.98 

0.77 

6.99 

3 

0.97 

0.76 

7.44 

Unnamed 

1 

0.91 

0.63 

4.30 

2 

0.92 

0.67 

5.47 

3 

1  .01 

0.74 

7.64 

Shriner 

1 

0.86 

0.72 

6.32 

2 

0.92 

0.74 

7.86 

3 

0.96 

0.76 

6.46 

Crystal 

1 

0.96 

0.90 

6.78 

2 

0.94 

0.89 

6.99 

3 

0.94 

0.86 

8.09 

Allen 

1 

0.97 

0.80 

6.91 

2 

1  .03 

0.73 

7.32 

3 

1.10 

0.32 

7  .  37 

Luii i  se 

1 

0.38 

0.61 

4  .  36 

Shadow 

1 

1.00 

0.88 

6.95 

George 

1 

0.96 

0.86 

7.02 

Marsh 

2 

1.09 

0.73 

5.20 

Tipsoo 

2 

0.91 

0  .  85 

7.34 

3 

0.96 

0.90 

7.37 

Frozen 

2 

1  .26 

0.78 

4.77 

Eunice 

2 

0.90 

0.66 

4.02 

Reflection 

0.95 

0.74 

7.72 

Eleanor 

1.04 

0.74 

7.28 

Fan 

3 

1  .09 

0.61 

3.32 

max 

1.26 

0.90 

8.09 

min 

0.86 

0.61 

3.32 

avg 

0.98 

0.76 

6.23 

std 

0.08 

0.09 

1.33 

Sea  Salt  Corr.  Major  Ion  Cone.,  ueq/1 


2  + 

2+ 

- 

- 

p  - 

- 

Lake     Tr 

ip 

Ca 

Mg 

Xa 

K 

SO 

4 

HC0 

o 

o 

Golden 

1 

41 

84 

2 

79 

19 

07 

1 

37 

4 

09 

58 

09 

2 

44 

90 

o 

78 

20 

83 

1 

37 

6 

84 

58 

29 

3 

48 

20 

4 

40 

23 

68 

1 

48 

8 

38 

61 

36 

Mowich 

1 

71 

00 

4 

42 

14 

43 

0 

45 

4 

01 

79 

23 

2 

75 

33 

4 

28 

11 

91 

0 

71 

3 

74 

84 

60 

3 

73 

42 

4 

02 

13 

51 

0 

69 

4 

55 

81 

88 

un -named 

1 

39 

42 

3 

21 

28 

15 

1 

39 

6 

80 

56 

01 

2 

50 

31 

4 

38 

27 

35 

1 

05 

8 

43 

66 

64 

56 

75 

4 

61 

28 

66 

0 

67 

9 

44 

81 

i  ° 

Shriner 

1 

107 

46 

8 

38 

29 

99 

5 

57 

12 

13 

120 

73 

2 

124 

50 

9 

93 

21 

91 

2 

11 

17 

52 

132 

14 

3 

110 

10 

9 

14 

20 

32 

0 

96 

16 

50 

122 

29 

Crystal 

1 

96 

02 

4 

59 

19 

06 

3 

96 

2 

92 

117 

59 

2 

112 

35 

4 

71 

17 

30 

1 

39 

3 

65 

126 

30 

3 

112 

73 

4 

35 

13 

91 

1 

77 

3 

95 

121 

50 

Allen 

1 

384 

88 

25 

75 

32 

35 

0 

39 

64 

04 

386 

69 

2 

428 

14 

29 

12 

31 

93 

1 

99 

112 

94 

413 

on 

o 

O 

397 

97 

23 

02 

32 

28 

0 

70 

120 

70 

4  1  1 

66 

Louise 

1 

23 

09 

0 

98 

12 

61 

0 

51 

0 

00 

29 

33 

Shadow 

1 

32 

13 

3 

66 

12 

35 

1 

37 

2 

25 

19 

51 

George 

1 

166 

92 

9 

50 

27 

08 

2 

97 

1 

60 

198 

10 

Marsh 

2 

32 

87 

4 

15 

20 

27 

2 

64 

14 

04 

53 

77 

Tipsoo 

2 

232 

90 

7 

63 

23 

90 

o 

21 

10 

31 

234 

77 

3 

207 

68 

7 

42 

26 

40 

2 

25 

9 

12 

228 

57 

Frozen 

2 

10 

12 

0 

69 

10 

11 

0 

54 

9 

98 

19 

o  ■> 

yj   J. 

Eunice 

2 

33 

98 

1 

28 

12 

75 

0 

47 

1 

56 

38 

58 

Reflection 

41 

59 

4 

92 

18 

82 

0 

66 

6 

66 

58 

08 

Eleanor 

3 

113 

31 

13 

71 

oo 

98 

0 

48 

34 

87 

1  oo 

J  O*  > 

86 

Fan 

o 

O 

21 

67 

1 

90 

15 

60 

1 

79 

12 

28 

O  O 

65 

max 

428 

14 

29 

12 

32 

35 

5 

57 

120 

70 

413 

00 

min 

10 

12 

0 

69 

10 

11 

0 

45 

0 

00 

19 

31 

avg 

113 

50 

7 

44 

21 

06 

1 

53 

17 

72 

126 

26 

std 

111 

97 

7 

44 

6 

74 

1 

13 

29 

59 

109 

05 

79 


Ratios 

HC03/ 

HC03/ 

'hco~~ 

S04, 

S04  m 

<§"• 

S04 

Tot.S04 

Lake 

Trip 

Cat. 

S04 

Cond. 

Anal . 

Bkg. 

Sea+Bkg 

/Sea 

Golden 

1 

0.89 

14.21 

3.87 

3.66 

24 

46 

-16 

71 

-20.80 

2 

0.82 

8.53 

5.55 

3.62 

26 

58 

-16 

13 

-22.96 

3 

0.79 

7.33 

5.55 

4.24 

29 

89 

-17 

27 

-25 . 65 

Mowich 

1 

0.88 

19.77 

5.32 

3.24 

21 

67 

-14 

43 

-18.43 

2 

0.92 

22.63 

6.99 

4.32 

19 

78 

-11 

72 

-15.46 

3 

0.89 

18.00 

7.44 

3.86 

21 

17 

-12 

76 

-17.31 

Unnamed 

1 

0.78 

8.24 

4.34 

3.65 

33 

38 

-22 

93 

^29.73 

2 

0.80 

7.91 

5.51 

3.42 

33 

-21 

40 

-29.91 

3 

0.89 

8.60 

7  .  65 

4.06 

34 

35 

-21 

36 

30.80 

Shriner 

1 

0.80 

9.95 

6.32 

2.  13 

45 

99 

-31 

73 

-43. 8G 

2 

0.83 

7.54 

7.87 

2.62 

36 

12 

-15 

98 

3 

0.87 

7.41 

6.47 

3.28 

32 

01 

-12 

23 

-23.73 

Crystal 

1 

0.95 

40.25 

6.80 

0.81 

32 

07 

-28 

35 

-31 .27 

2 

0.93 

34.57 

7.00 

0.55 

29 

23 

-25 

02 

-28.68 

3 

0.92 

30.74 

8.10 

1  .34 

26 

21 

-20 

92 

-24.87 

Allen 

1 

0.37 

6.04 

6.91 

4.26 

70 

19 

-1 

90 

-65.94 

2 

0.85 

3.70 

7  .  32 

4.23 

75 

08 

12 

08 

-70.35 

O 

0.90 

3.41 

7.36 

4.44 

71 

32 

53 

-  66.33 

Louise 

1 

0.79 

N/A 

4  .  93 

2.08 

15 

29 

-13 

21 

-13.21 

Shadow 

1 

1  .00 

22.00 

6.97 

0.79 

16 

95 

-13 

91 

-16.16 

George 

1 

0.96 

123.62 

7.02 

4.87 

45 

92 

-39 

45 

-41 .05 

Marsh 

2 

0.90 

3.33 

5.27 

2.41 

26 

24 

-9 

79 

-23.33 

Tipsoo 

2 

0.88 

21  .71 

7.34 

1  .20 

47 

84 

-35 

82 

-46.64 

3 

0.94 

25.06 

7.37 

1.35 

48 

01 

-37 

54 

-46.66 

Frozen 

2 

0.90 

1  .94 

5.28 

0.08 

11 

62 

-1 

57 

-11 .55 

Eunice 

2 

0.80 

24.77 

4.06 

2.61 

16 

39 

-12 

23 

-13.79 

Reflect. 

ion    3 

0.88 

8.72 

7.74 

2.11 

23 

67 

-14 

91 

-21 .56 

Eleanor 

3 

0.88 

3  .  84 

7.28 

3.96 

35 

89 

2 

94 

-31 .93 

Fan 

0.82 

2.74 

3.38 

2.50 

19 

51 

-  J 

/  3 

- 17 .01 

max 

1.00 

123.62 

8.  10 

4.87 

75 

08 

53 

82 

-11 .55 

min 

0.78 

0.00 

3.58 

0.08 

11 

62 

-39 

45 

-70. 85 

avg 

0.87 

17.14 

6.32 

2.82 

33 

47 

-12 

94 

-30.66 

std 

0.06 

22.64 

1.28 

1  .33 

16 

26 

19 

48 

15.80 

3C 


Table   A3-2: 


Aluminum,    ug/1 


PH 


Trip 

Total 

Organic 

Lab 

Lake 

Field 

Lab 

ile 

Golden 

1 

6.77 

7.75 

90.00 

NA 

NA 

2 

6.20 

7.40 

88.72 

NA 

NA 

3 

6.28 

7.60 

89.00 

NA 

NA 

Mowich 

1 

5.93 

7.04 

2.66 

NA 

NA 

2 

5.56 

7.09 

0.79 

NA 

NA 

3 

6.68 

6.92 

2.01 

NA 

NA 

Unnamed 

1 

5.80 

6.47 

1  .14 

NA 

NA 

2 

5.78 

6.35 

2.39 

NA 

NA 

3 

6  .  35 

6.75 

4.51 

NA 

NA 

Shriner 

1 

6.  17 

7.  16 

14.00 

NA 

NA 

2 

G.09 

7.  18 

12.45 

NA 

NA 

3 

5.73 

6.78 

14.47 

NA 

NA 

Crystal 

I 

5.71 

7.25 

1  .08 

NA 

NA 

2 

6.76 

6.86 

4.01 

NA 

NA 

3 

6.72 

6.82 

10.70 

NA 

NA 

Allen 

1 

6.64 

7.61 

2.37 

NA 

NA 

2 

7.34 

7.68 

3.01 

NA 

NA 

3 

6.66 

7.55 

4.31 

NA 

NA 

Louise 

1 

6.10 

6  .  03 

2.45 

NA 

NA 

Shadow 

1 

6.50 

6.61 

2.81 

NA 

NA 

George 

1 

6.30 

7.25 

0.89 

NA 

NA 

Marsh 

2 

7.31 

5.98 

3.29 

NA 

NA 

Tipsoo 

2 

7.20 

7.07 

78.47 

NA 

NA 

Frozen 

2 

5.55 

5.70 

0.33 

NA 

NA 

Eunice 

2 

5 .  53 

6.30 

1  .83 

NA 

NA 

Reflectio   3 

6.  14 

6.57 

6.26 

NA 

NA 

Eleanor 

3 

7.06 

7.09 

2.19 

NA 

NA 

max 

7  .  34 

7.75 

90.00 

NA 

NA 

min 

5.53 

5.70 

0.79 

NA 

NA 

avg 

6.33 

6.92 

16.56 

NA 

NA 

std 

0.53 

0.53 

29.49 

NA 

NA 

Surface 

+  Depth 

max 

7.43 

7.75 

90.00 

17 

.90 

70 

.26 

min 

0.00 

0.00 

0.00 

0 

.00 

0 

.00 

avg 

5.78 

6.50 

15.03 

2 

.99 

11 

.95 

std 

1.88 

1.52 

25.43 

4 

.39 

19 

.46 

ANC 


Anions,  mg/1 


Lake 


Golden 


Mowich 


Unnamed 


Shriner 


Crystal 


Allen 


Trip 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 
o 

1 

2 
3 
1 
2 


Louise 

Shadow 

George 

Marsh 

Tipsoo 

Frozen 

Eunice 

Reflection 

Eleanor 

max 
m  i  n 
avg 
std 


max 
min 
avg 
std 


ueq/1 


680.00 

673.20 

693.50 

81  .69 

93.65 

90.10 

99.60 

92.20 

86.98 

120.68 

131 .38 

120.70 

134.60 

122.30 

132.50 

438.42 

416.70 

400.94 

35.40 

45  .  80 

213.00 

58.95 

235.00 

10.30 

35.46 

62.00 

138.60 


mg/1 

as 

CaCO 

34.00 
33  .  66 
34.68 


08 
68 
51 
98 
61 
35 


6.03 

6.57 

6.04 

6.73 

6.  12 

6.63 

21.92 

20.84 

20  .  05 

1  .77 

2.29 

10.65 

2.95 

11.75 

0.52 

1.77 

3.10 

6.93 


693.50  34.68 

10.30  0.52 

201.62  10.08 

202.21  10.11 


693 . 50  34 . 68 

0.00  0.00 

155.03  7.63 

164.11  8.20 


2- 

- 

- 

2- 

- 

co 

HC0 

CI 

SO 

F(ug/ 

3 

3 

4 

0 

05749 

41 

48 

0 

51 

0 

90 

37 

99 

0 

02543 

41 

07 

0 

89 

2 

83 

26 

80 

0 

04151 

42 

30 

1 

12 

1 

06 

24 

52 

0 

00135 

4 

99 

0 

59 

.  0 

31 

11 

20 

0 

00173 

5 

72 

0 

55 

0 

46 

10 

04 

0 

00113 

5 

50 

0 

51 

0 

49 

10 

39 

0 

00044 

6 

10 

0 

57 

0 

31 

15 

75 

0 

00031 

5 

65 

0 

46 

0 

46 

14 

33 

0 

00074 

5 

32 

0 

53 

0 

65 

11 

77 

0 

00262 

7 

37 

0 

27 

0 

68 

15 

75 

0 

00299 

8 

02 

0 

24 

96 

1 5 

39 

0 

00109 

7 

37 

0 

4  2 

0 

91 

1  ° 

J  KJ 

87 

0 

00360 

8 

21 

1 

07 

0 

20 

L  O 

27 

0 

00133 

7 

47 

0 

10 

0 

22 

15 

88 

0 

00132 

8 

09 

0 

22 

0 

26 

13 

69 

0 

02685 

26 

74 

0 

51 

4 

06 

20 

42 

0 

02999 

25 

42 

0 

59 

6 

03 

17 

93 

0 

02139 

24 

46 

0 

56 

6 

17 

20 

fin 

0 

00006 

2 

22 

0 

41 

0 

20 

10 

29 

0 

00028 

2 

31 

0 

37 

0 

30 

15 

75 

0 

00570 

13 

00 

0 

62 

0 

85 

13 

35 

0 

00009 

3 

66 

0 

23 

0 

87 

20 

25 

0 

00415 

14 

34 

0 

15 

0 

59 

24 

98 

0 

00001 

0 

75 

0 

05 

0 

29 

9 

96 

0 

00011 

2 

19 

0 

66 

0 

20 

6 

79 

0 

00035 

3 

80 

0 

27 

0 

43 

11 

20 

0 

00257 

8 

46 

0 

61 

1 

87 

19 

88 

0 

05749 

42 

30 

1 

12 

6 

17 

37 

99 

0 

00001 

0 

75 

0 

05 

0 

20 

6 

79 

0 

00869 

12 

31 

0 

49 

I 

21 

16 

60 

0 

01467 

12 

32 

0 

25 

1 

62 

6 

48 

0 

05749 

42 

30 

1 

12 

6 

17 

37 

99 

0 

00000 

0 

00 

0 

00 

0 

00 

0 

00 

0 

00574 

9 

33 

0 

40 

1 

11 

14 

54 

0 

01121 

9 

99 

0 

23 

1 

60 

6 

27 

Cations,  mg/1 


2  + 

2+ 

+ 

- 

Lake 

Trip 

Ca 

Mg 

Na 

K 

SiO 

Golden 

1 

10.72 

1.17 

2.24 

0.58 

18 

57 

2 

10.43 

1.09 

2.04 

0.71 

19 

11 

3 

11.  10 

1.13 

1.95 

0.68 

19 

11 

Mowich 

1 

1.56 

0.19 

0.51 

0.03 

2 

96 

2 

1.52 

0.  18 

0.48 

0.06 

3 

09 

3 

1.68 

0.19 

0.50 

0.03 

3 

30 

Unnamed 

1 

1.20 

0.23 

0.85 

0.  10 

7 

73 

2 

1.32 

0.23 

0.93 

0.16 

6 

65 

O 

1  .  11 

0.20 

0.85 

0.  12 

5 

98 

Shriner 

1 

2.07 

0.26 

0.79 

0.16 

3 

70 

2 

2.49 

0.33 

0.62 

0.08 

o 
o 

16 

3 

2.22 

0.32 

0.62 

0  .  05 

2 

25 

Crystal 

1 

2.26 

0.  16 

0.48 

0.  13 

66 

2 

2.30 

0.  15 

0.43 

0.06 

3 

55 

o 

2.27 

0.  15 

0.38 

0.05 

3 

29 

Allen 

1 

9.22 

0.98 

0.97 

0  .  05 

6 

25 

2 

8.59 

0.93 

0.93 

0.08 

3 

32 

3 

8.35 

0  .  83 

0.96 

0  .  05 

3 

08 

Louise 

1 

0.60 

0.07 

0.47 

0.04 

O 

6o 

Shadow 

1 

0.64 

0.  13 

0.32 

0  .  08 

74 

George 

1 

3.86 

0.35 

0.86 

0.  13 

4 

24 

Marsh 

2 

1  .00 

0.  19 

0  .  58 

0.  10 

4 

20 

Tipsoo 

2 

4.88 

0.23 

0.60 

0.  12 

2 

33 

Frozen 

2 

0.17 

0.03 

0.24 

0.04 

O 

36 

Eunice 

2 

0.75 

0.07 

0.44 

0.03 

1 

64 

Reflect. 

ion 

0.85 

0.19 

0.55 

0.04 

3 

80 

Eleanor 

3 

2.41 

0.48 

0.73 

0.02 

4 

69 

max 

11  .  16 

1  .17 

2  24 

0.71 

19 

11 

in  in 

0.  17 

0.03 

0.24 

0.02 

1 

64 

avg 

3.54 

0.39 

0.79 

0.  14 

5 

53 

std 

3.49 

0.36 

0.50 

0.  19 

4 

93 

max 

11  .  16 

1  .  17 

2.24 

0.71 

19 

11 

min 

0.00 

0.00 

0.00 

0.00 

0 

00 

avg 

2.78 

0.30 

0.66 

0.  10 

4 

29 

std 

2.99 

0.30 

0.39 

0.  13 

3 

61 

Calc  . 

Total 

Total 

Ra  t  i  o 

DOC 

COIK 

I. 

Ionic 

Cations 

Anions 

Tot .Cat 

Lake 

Trip 

mg/1 

umho/cm 

Str. 

ueq/1 

ueq/1 

/An . 

mole/1 

Golden 

1 

20.63 

83 

00 

0 

00106 

734.01 

715.02 

1  .03 

2 

22.31 

77 

50 

0 

00107 

708.32 

758.69 

0.93 

3 

24.77 

75 

05 

0 

00109 

741  .76 

748.49 

0.99 

Mowich 

1 

1  .20 

16 

50 

0 

00016. 

115.07 

105.49 

1  .09 

2 

0.50 

11 

90 

0 

00017 

112.13 

119.14 

0.94 

3 

1.20 

12 

50 

0 

00017 

120.68 

115.32 

1  .05 

Unnamed 

1 

0.25 

15 

20 

0 

00016 

116.79 

122.94 

0.95 

2 

0.30 

13 

00 

0 

00017 

127.82 

115.47 

1.  11 

3 

0.28 

10 

60 

0 

00016 

111 .27 

115.98 

0.96 

Shriner 

1 

0.97 

18 

90 

0 

00022 

161.10 

143.36 

1  .  12 

2 

1  .30 

16 

80 

0 

00026 

177.32 

159.03 

1  .  12 

3 

1  .50 

19 

10 

0 

00024 

162.81 

152.37 

1  .07 

Crystal 

1 

0.90 

19 

00 

0 

00022 

147.69 

169.65 

0.87 

2 

1.20 

18 

70 

0 

00020 

145.20 

130.59 

1.11 

3 

1.30 

13 

50 

0 

00021 

141 .21 

144.85 

0.97 

Allen 

1 

0.73 

65 

10 

0 

00087 

576.02 

538.53 

1  .07 

o 

0.75 

57 

10 

0 

00087 

540.68 

559 . 64 

0.97 

3 

0.75 

55 

00 

0 

00085 

524.72 

54  6.36 

0.96 

Louise 

1 

0.45 

7 

50 

0 

00007 

56.53 

51  .70 

1  .09 

Shadow 

1 

1.50 

11 

00 

0 

00009 

57  .  93 

63  .  34 

0.91 

George 

1 

0.83 

29 

00 

0 

00038 

258.59 

248.83 

1  .04 

Marsh 

2 

0.26 

10 

40 

0 

00013 

92.64 

86.04 

1  .08 

Tipsoo 

2 

3.00 

32 

00 

0 

00041 

287.37 

252.87 

1.14 

Frozen 

0.23 

3 

47 

0 

00003 

21  .87 

18.30 

1  .20 

Eunice 

2 

0.01 

9 

20 

0 

00008 

62.24 

58  .  60 

1  .06 

Reflect. 

loii 

2.20 

7 

35 

0 

00011 

81  .72 

79.  14 

1  .03 

Eleanor 

3 

0.67 

19 

20 

0 

00029 

189.77 

195.72 

0.97 

max 

24.77 

33 

00 

0 

00109 

741 .76 

758.69 

1  .  20 

min 

0.01 

47 

0 

00003 

21  .37 

18.30 

0.87 

avg 

3.33 

26 

95 

0 

00036 

243.47 

241.31 

1.03 

std 

6.85 

23 

64 

0 

00034 

222.02 

226.18 

0.08 

Surface 

+  Depth 

max 

24.77 

83 

00 

0 

00109 

741  .76 

758.69 

1.20 

min 

0.00 

0 

00 

0 

00000 

0.00 

0.00 

0.00 

avg 

2.17 

22 

08 

0 

00028 

191 .58 

189.77 

0.96 

std 

4.69 

19 

93 

0 

00029 

186.82 

191.13 

0.26 

Ratio    Ratio   Ratio 
An/Cat  ANC/Cat  ANC/Cond 
Lake     Trip 


Golden 

1 

0.97 

0.93 

8.  19 

2 

1.07 

0.95 

8.69 

3 

1.01 

0.93 

9.24 

Mowich 

1 

0.92 

0.71 

4.95 

2 

1  .06 

0.84 

7.87 

3 

0.96 

0.75 

7.21 

Unnamed 

1 

1.05 

0.85 

6.55 

2 

0.90 

0.72 

7.09 

3 

1  .04 

0.78 

8.2] 

Shriner 

1 

0.89 

0.75 

6.38 

0.89 

0.74 

7.32 

3 

0.94 

0.74 

6.32 

Crystal 

1 

1.15 

0.91 

7.08 

2 

0.90 

0.84 

6.54 

3 

1  .03 

0.94 

9.81 

Allen 

1 

0.93 

0.76 

6.73 

2 

1.04 

0.77 

7.30 

3 

1  .04 

0.76 

7.29 

Lou a se 

i 

0.91 

0  .  63 

4  .72 

Shadow 

1 

1  .09 

0.79 

4.  16 

George 

1 

0.96 

0.32 

7.34 

Marsh 

2 

0.93 

0.64 

5.67 

Tipsoo 

2 

0.88 

0.82 

7.34 

Frozen 

2 

0.84 

0.47 

2.97 

Eunice 

2 

0.94 

0.57 

3.85 

Reflection 

0.97 

0.76 

8.44 

Eleanor 

1  .03 

0.73 

7.22 

max 

1.  15 

0  .  95 

9.81 

min 

0.84 

0.47 

2.97 

avg 

0.98 

0.77 

6.85 

std 

0.08 

0.11 

1  .60 

Surface  + 

Depth 

max 

1  .26 

0  .  95 

9.81 

min 

0.00 

0.00 

0.00 

avg 

0.92 

0.72 

6.13 

std 

0.25 

0.21 

2.  12 

Sea  Salt  Corr.  Major  Ion  Cone.,  ueq/1 


2  + 

2  + 

- 

- 

c  - 

- 

Lake     Tr 

ip 

Ca 

Mg 

Na 

K 

SO 

HC0 

4 

3 

Golden 

1 

534 

11 

40 

59 

89 

61 

14 

61 

14 

61 

679 

90 

2 

519 

40 

36 

51 

74 

65 

17 

62 

51 

96 

673 

12 

3 

555 

23 

36 

76 

67 

14 

16 

66 

13 

08 

693 

40 

Mowich 

1. 

77 

01 

4 

68 

12 

88 

0 

50 

1 

87 

81 

77 

2 

75 

27 

4 

57 

12 

38 

1 

18 

5 

21 

93 

71 

3 

83 

06 

5 

07 

13 

83 

0 

37 

6 

15 

90 

20 

Unnamed 

1 

58 

91 

6 

05 

28 

00 

2 

27 

1 

85 

99 

92 

2 

65 

14 

6 

50 

33 

42 

3 

73 

5 

96 

92 

63 

3 

54 

79 

5 

43 

28 

63 

2 

77 

9 

38 

87 

14 

Shriner 

1 

102 

93 

8 

4  3 

29 

91 

3 

88 

12 

00 

120 

72 

2 

123 

99 

10 

96 

°2 

96 

1 

96 

18 

08 

131 

■12 

3 

110 

02 

9 

94 

20 

32 

1 

07 

15 

66 

120 

84 

Crystal 

1 

111 

19 

1 

70 

o 
O 

89 

2 

75 

0 

00 

134 

63 

2 

114 

84 

4 

97 

17 

06 

1 

44 

3 

78 

122 

42 

3 

112 

99 

4 

52 

13 

00 

1 

23 

3 

69 

132 

63 

Allen 

1 

459 

49 

33 

64 

34 

04 

1 

09 

80 

47 

438 

37 

2 

427 

96 

31 

73 

31 

44 

1 

78 

120 

84 

416 

65 

3 

416 

10 

29 

87 

32 

79 

0 

94 

124 

03 

400 

90 

Louise 

1 

29 

61 

0 

89 

13 

32 

0 

85 

0 

92 

36 

32 

Shadow 

1 

31 

52 

26 

3 

10 

1 

78 

o 
o 

46 

0-1 

George 

1 

191 

53 

10 

49 

27 

58 

2 

96 

12 

83 

213 

02 

Marsh 

2 

49 

65 

5 

77 

21 

03 

2 

46 

15 

90 

59 

99 

Tipsoo 

2 

243 

14 

7 

87 

23 

68 

3 

05 

11 

03 

235 

04 

Frozen 

2 

3 

65 

0 

74 

9 

43 

0 

94 

5 

64 

12 

29 

Eunice 

2 

36 

60 

0 

00 

8 

83 

0 

29 

0 

00 

35 

96 

Reflection 

41 

87 

5 

82 

19 

52 

0 

89 

6 

83 

62 

26 

Eleanor 

3 

119 

23 

15 

10 

22 

31 

0 

21 

34 

00 

138 

66 

max 

555 

23 

40 

59 

89 

61 

17 

62 

124 

03 

693 

40 

min 

8 

65 

0 

00 

3 

09 

0 

21 

0 

00 

12 

29 

avg 

176 

08 

12 

29 

26 

68 

3 

31 

21 

45 

201 

85 

std 

173 

92 

12 

59 

19 

93 

4 

72 

33 

25 

202 

00 

Surface  +  Depth 

max           555.23  40.59  89.61  17.62  124.03  693.40 

min             0.00  0.00  0.00  0.00  0.00  0.00 

avg           136.09  9.55  22.60  2.30  20.01  155.31 

std           148.67  10.52  15.07  3.38  32.68  163.90 


Ratios 

HC03/ 

HC03/ 

HC0  / 

S04, 

S04  Bkg. 

S04 

Tot .S04 

Lake 

Trip 

Cat. 

S04 

Cond. 

Anal  . 

Bkg. 

Sea+Bkg 

/Sea 

Golden 

1 

1.00 

46.55 

8.19 

4.03 

155 

95 

-137.31 

-151 

92 

2 

1  .04 

12.95 

8.69 

7.03 

142 

30 

-83.32 

-135 

28 

3 

1.03 

53.00 

9.24 

8.88 

137 

08 

-115.11 

-128 

19 

Mowich 

1 

0.86 

43.72 

4.96 

4.67 

20 

73 

-14.20 

-16 

07 

2 

1  .00 

18.00 

7.88 

4.32 

20 

74 

-11 .21 

-16 

42 

3 

0.88 

14.67 

7.22 

4.05 

22 

13 

-11 .93 

-18 

08 

Unnamed 

1 

1.05 

53.90 

6.57 

4.52 

36 

12 

-29.75 

-31 

60 

2 

0.85 

15.55 

7.  13 

3.62 

43 

60 

-34.02 

-39 

98 

o 
O 

0.95 

9.29 

8.22 

4.  16 

36 

32 

-23.28 

-32 

66 

Shriner 

1 

0.83 

10.06 

6.39 

2.  16 

43 

81 

-29.65 

-41 

66 

2 

0.82 

7.27 

7.82 

1  .91 

37 

06 

-17.07 

-35 

15 

3 

0.85 

7.71 

6.33 

3.28 

32 

19 

-13.24 

-28 

90 

Crystal 

1 

1  .  13 

N/A 

7.09 

4.16 

16 

80 

-12.64 

-12 

64 

2 

0.89 

32.43 

6.55 

0.81 

29 

28 

-24.70 

-28 

47 

3 

1.01 

35.98 

9.82 

1.74 

24 

81 

-19.39 

-23 

08 

Allen 

1 

0.83 

5.45 

6.73 

4.06 

79 

52 

5.01 

-75 

46 

2 

0.85 

3.45 

7.30 

4.63 

74 

59 

50  .  87 

-69 

97 

3 

0.84 

3.23 

7.29 

4.44 

73 

86 

54.61 

-69 

42 

Louise 

1 

0.80 

39  .  56 

1 .  34 

3.25 

17 

42 

-13.26 

-14 

13 

Shadow 

1 

1.03 

13.88 

4.  19 

2.93 

13 

01 

-6.76 

-10 

08 

George 

1 

0.92 

16.61 

7.35 

4.87 

48 

73 

-31 .03 

-43 

86 

Marsh 

2 

0.76 

3.77 

5.77 

2.21 

28 

47 

-10.36 

-26 

26 

Tipsoo 

2 

0.85 

21  .30 

7.35 

1  .21 

49 

32 

-37 . 08 

-48 

11 

Frozen 

2 

0.62 

2.  18 

3.54 

0.40 

11 

21 

-5.  18 

-10 

81 

Eunice 

2 

0.79 

N/A 

3.91 

4.16 

12 

41 

-8.25 

-8 

25 

Reflection 

0.91 

9.  11 

8.47 

2.  12 

24 

70 

-15.74 

-22 

58 

Eleanor 

3 

0.88 

4.08 

7.22 

4.83 

34 

61 

4.22 

-29 

78 

max 

1  .13 

53.90 

9.82 

8.88 

155 

95 

54  .61 

-8 

25 

min 

0.62 

0.00 

3  .  54 

0.40 

11 

21 

-137.31 

-  1  5  1 

92 

avg 

0.90 

17.91 

6.89 

3.65 

46 

94 

-21  .84 

-43 

29 

std 

0.11 

16.53 

1.54 

1.78 

39 

16 

38  .  65 

38 

13 

Surface 

+  Depth 

max 

1.13 

123.62 

9.82 

8.88 

155 

95 

54.61 

15 

80 

min 

0.00 

0.00 

0.00 

0.00 

0 

00 

-137.31 

-151 

92 

avg 

0.83 

18.45 

6.17 

3.00 

38 

92 

-15.23 

-34 

76 

std 

0.23 

23.56 

2.09 

1.74 

30 

24 

31.03 

30 

00 

Appendix 


Nomenclature 


[i] 


Partial  pressure  of  atmospheric  C0? 
Molar  concentration  of  any  species,  i 


Z  [cations]  :  equivalent  sum  of  major  cations 


D.O. 
PH 


=  (2[Ca2~]+2[Mg2^]-[Na~]-[K~] 


Dissolved  oxygen 
Negative  log  [H  ] 


[ANC]___ , carbonate  ;   Acid  neutralizing  capacity;   equivalent? 

of  acid  required  to  titrate  a  water  to  the  El«C0o  equivalence 
point  (f=0)  for  the  carbonate  system;  equal  to  ^Alkalinity] 
for  carbonate  buffered  waters 


where  [Alkalinity]  =  [HC0o  ]-2[C0o-  3+[0H~]-[K  ] 
=  [HC0o  ]  for  pH  range  5.5  to  "   5 

:    First,  acid  dissociation  constant  .  carbonate  system 

:    Henry's  Law  constant  for  C0o(g)  solubility  in  water 

:     Composite  first  acid  dissociation  constant, 

carbonate  system  =  K  K,T 

a!  h 

:    Acid  dissociation  constant  for  mode]  organic  acid 
Solubility  constant  for  solid  phase  dissocial 
Negative  log  of  solubility  constant 
Change  in  (loss  of) 

Correlation  coefficient  squared  (coefficient  of  determination 

Empirical  constant  equaling  =  log  b'-log  K    log  orn 

a1     '  ^>0 

Empirical  constants 


K 
a 

or 

K 

X  I ) 

so 

A 

2 

r 

a 

b. 

b' 

